Renormalization group for 2D fermions with a flat Fermi surface 
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We present a renormalization group analysis of two-dimensional interacting fermion systems with 
a closed and partially flat Fermi surface. Numerical solutions of the one-loop flow equations show 
that for a bare local repulsion, the system evolves through three different regimes as the typical 
energy is lowered : a high-energy transient with a strong competition between particle-particle and 
particle-hole channels, an intermediate regime with dominant spin density wave correlations, and 
finally a "hot spot" regime exhibiting d-wave superconductivity. We study, mostly by analytical 
methods, how this flow pattern depends on the number A*' of Fermi surface patches used in the 
numerical solution. This clearly indicates that the final regime requires vanishingly small microscopic 
interactions, for the one-loop approximation to be valid, as A*' is going to infinity. 
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I. INTRODUCTION 

Over the past years, several groups have been 
using renormalization group (RG) methods to 
study |*i¥jft"|4HWi>W'^^l (2D) interacting electron 
systems Jd'ElijQijLllirlalj One of the main advantage of 
these approaches is that in principle they are not biased 
towards any special type of instability. Charge and 
spin ordering and superconductivity are treated on the 
same footing, which is highly desirable in the modelling 
of many materials such as oxide superconductors or 
organic conductors. The quoted works often deal with 
rather complex situations, since they are inspired by 
real systems. In fact, this is another advantage with 
the RG schemes based on successive elimination of 
high-energy modes : basically no assumption (such as 
renormalizability) has to be made on the underlying 
microscopic model. It has therefore been possible 
to stucb/pthe effect of Van-Hoiic ainftularities, in the 
absenccOO'El or in the presencoJO lIulIIj'Ej of partial or 
perfect Fermi surface nesting. Some very encouraging 
results have been obtained in this way, showing how 
d-wave superconductivity iflMty occur in a doped Mott 
insulating antiferromagnetEicrB or even how the metallic 
character of such systems gradually disappears in 
momentum space, starting from saddle points to extend 
progressively over the ftjbole Fermi line, as the typical 
energy scale is reducedrllj 

The present work is motivated by the desire to analyze 
in detail a slightly simpler situation, which still captures 
most of the complexity of some two-dimensional electron 
systems. We shall consider the case of a closed Fermi 
line which exhibits some flat segments. In this situation, 
we have perfect Fermi surface nesting, but no singularity 
in the single-particle density of states. Therefore, quan- 
tum corrections to the effective interaction produce the 



same logarithmic divergence in the particle-particle (or 
Cooper) channel as in the particle-hole (or Peierls) chan- 
nel. As recognized already long ago by the Soviet School, 
this leads to a rather intricate competition between var- 
ious possible instabilities. The corresponding "parquet" 
diagrams have been investigate d . in , various situations : a 
semi-metal in a magnetic field ,F.HH the itinerant antifer- 
romagnetic state of chromium,Ej or -the phase diagram 
of quasi one-dimensional conductors. La These important 
works have discussed the existence of two types of so- 
lutions, referred to as "fixed poles" and "moving poles" 
according to whether all the effective couplings diverge 
at a single energy scale or instead decouple into various 
classes which each have their own instability scale. The 
moving poles occur mostly in situations where a single 
channel emerges in the low-energy regime, after a com- 
plicated transient at higher energies. This corresponds 
to the simple picture of particle-hole or particle-particle 
bound state formation. Since the total momentum of 
these bound states is a good quantum number in a pure 
system, we do get decoupled poles, one for each value 
of this momentum. For a nearly square Fermi surface, 
and with repulsive locaLinteractions, this scenario has 
already been supportedJi3 and confirmed more rpeently 
thanks to more extensive numerical calculations. ta The 
mobile poles then occur in the particle-hole channel, with 
the formation of a spin density wave. 

We have recently analyzed numerically the same prob- 
lem, with special focus on the relative evolution of variopi 
couplings in the immediate vicinity of the singularity.E3 
This was done using differential equations for the nor- 
malized couplings (i.e. the othogonal projection of the 
coupling vector onto the unit sphere in coupling space). 
To our surprise, we found that the spin density wave 
regime is not the final stage of the RG flow (according 
to the one- loop approximation). Seen in terms of nor- 



malized couplings, this regime corresponds to an inter- 
mediate fixed point which is ultimately unstable towards 
a d-wave superconductor ! It was then very tempting 
to relate this remarkable behavior of the RG flow to the 
temperature evolution of physical properties of^upercon- 
ducting cuprates in the underdoped regime,llZI in which 
strong but local antiferromagnetic correlations built up 
in the pseudo gap state above the superconducting tran- 
sition. The present article aims at providing a detailed 
explanation for the occurence of this unforeseen RG flow. 
Unfortunately, we shall establish that its experimental 
relevance is severely limited by the fact it requires van- 
ishingly small couplings if the one-loop approximation we 
are using remains under good control. Numerical solu- 
tions of RG equations naturally introduce a discrete set 
of N patches instead of the continuous Fermi line. At 
finite N, we may observe the final fixed point provided 
the initial couplings are chosen in a finite interval whose 
width shrinks to zero as N goes to infinity. 





FIG. 1: (a) : The square with rounded corners shaped Fermi 
surface considered in this work. The distance between the 
origin and any flat region is denoted by K, and 20 is the 
angular size of the segments seen from the origin. We have < 
< 7r/4. (b) : Particle-hole pairs with vanishing excitation 
energy and a common total momentum Q. 



II. ONE-LOOP RENORMALIZATION GROUP 

ANALYSIS OF TWO-DIMENSIONAL MODELS 

WITH A NESTED FERMI SURFACE 



We should note that our work has some common fea- 
tures with a serJes of studies on coupled one-dimensional 
chains ,LSL3tjO and we shall borrow some of the tools dc- 
velnned in these articles, specially in the works by Lin et 
al.EJ'Ej However these investigations do not make the as- 
sumption of flat Fermi surfaces, so we get quite different 
results, simply because the allowed low-energy couplings 
are different. 



This paper is organized as follows. Sec. O describes the 
two-dimensional model, the one-loop RG equations and 
our numeripd results, in an expanded version of our pre- 
vious work.Ej In Sec. y, we analyze a simplified model of 
couplfijd chains (the same as in the work by Zheleznyak 
et alEj), focusing on generalized Luttinger liquid types 
of fixed_|points which were first mentioned in one of our 
works.t3 Sec. |^ discusses the fate of "moving poles" solu- 
tions as N is large but finite, so that perturbations from 
the Cooper channel are still able to affect the leading 
Peierls channel. In Sec. M we introduce special interme- 
diate fixed directions with large SU(7V) symmetry, and 
establish their physical properties, via the computation 
of some response functions. Sec. N^ discusses the local 
stability of these special directions in coupling space, and 
confirms the existence of unstable perturbations at fi- 
nite N, as observed numerically and expected on general 
grounds in Sec. ^A The ideas used in this work are then 
applied in Sec. VI]| to some other aspects of symmetry 
restoration under RG flows. An Appendix also contains 
some considerations regarding the interplay between RG 
and continuous global symmetries. 



A. CHOICE OF THE MODEL 

Let us consider a two-dimensional model of non- 
interacting electrons for which the Fermi surface exhibits 
the shape shown on Fig. 0(a), namely a square with 
rounded corners. The flat regions will be referred to 
as nested regions, since it is possible to excite a large 
family of particle-hole pairs with vanishing excitation en- 
ergy, provided the particle and the hole are created in 
the vicinity of the Fermi surface. A remarquable fea- 
ture of these Fermi surfaces is that they yield continuous 
one-parameter families of wave vectors (see Fig. |^(&)) for 
which the static charge or spin density response func- 
tions diverge in the zero temperature limit. With the 
notations of Fig. |l|(a), these special wave vectors Q are 
given by Q = {±2K,ky), or Q = {k^,±2K) with \k^\ 
and \ky\ freely chosen in the interval [0, 2_ft'tan8]. We 
couldn't flnd simple tight binding Hamiltonians which 
would lead to precisely this Fermi surface. But the tight 
binding Hamiltonian with nearest neighbor hopping on a 
square lattice is well known to generate a perfect square 
Fermi surface at half-filling. For a slightly less than 
half filled band, the Fermi surface is closed and remains 
highly anisotropic, with a diamond-like shape, although 
its curvature never vanishes. Experimentally, some very 
fiat Fermi surfaces are seen an. photoemission studies of 
Bi2Sr2CaCu20s compoundsHj Note that however some 
controversies do not seem to have been completely set- 
tled in the experimental community. AlterBajtive views 
are discussed for instance by Bogdanov et acJ and Feng 
et al£a 

Our simple model is expected to be valid as long as 
the temperature is not too small, so that the difference 
between the small curvature regions and the idealized flat 
segments lies in the thermally excited layer around the 
Fermi surface. 



Once the Fermi surface is given, we have to specify 
which interactions will be kept in an effective model, in 
order to capture the dominant processes, involving exci- 
tation energies already small with respect to the Fermi 
energy. Following the picture developed by Shankar,c3 
we consider all scattering processes with incoming mo- 
menta denoted by fci and fc2, and outgoing momenta by 
fca and k^. These are constrained by momentum con- 
servation : ki + k2 — ks + k4, and all the fc^'s have to 
lie on the Fermi surface. The allowed processes can be 
classified in four families. Three of these involve three 
continuous parameters each. By analogy toi-tlie standard 
notation for quasi one-dimensional systems pJ a first fam- 
ily which may be denoted by gi{ki,k2,k3,ki), involves 
all the fci's on the same flat segment of the Fermi sur- 
face (see Fig. 0(a)). Experience with one-dimensional 




FIG. 2: Examples of low-energy interaction processes, (a) : 
g4 process, (b) : g2 process, (c) : gi process, (d) and (e) : 
BCS interaction, (e) is peculiar because it involves the two 
pairs of flat FS segments. (/) shows one of the many ways to 
define "Right" and "Left" parts of the Fermi surface. 

systems suggest they won't play a leading role at low 
energies, since their effect is mostly to renormalize the 
single-particle Fermi velocity. We shall therefore discard 
them here, since we wish to focus on the possible static 
instabilities of our system. A second set of processes, 
denoted by 172(^1, fe, ^3, ^4), involves incoming particles 
on opposite parallel segments, and the transferred mo- 
mentum is parallel to the Fermi surface (see Fig. ||(&)). 
Note that for generic systems, namely anisotropic Fermi 
surfaces with non-vanishing curvature, the effective pa- 
rameters depend continuously on two variables. Indeed, 



in the low-energy limit, momentum conservation usually 
requires (^3,^4) = (^2,^1) or (fci,A:2). The presence of 
flat segments greatly enlarges the space of the available 
processes, and a third variable emerges, since the trans- 
ferred momentum may be varied at will along the flat di- 
rections. A third three parameter family generalizes the 
back-scattering processes of quasi one-dimensional sys- 
tem and will be denoted by 5i(fci, fc2, /ca, fc4). They are 
depicted on Fig. ||(c). Finally, we have the two parame- 
ter continuous family of BCS couplings for which the to- 
tal incoming momentum vanishes (see Fig. 0((i) and (e)). 
By contrast to the previous scattering configurations, the 
BCS couplings are able to connect one pair of flat regions 
to the other pair (Fig. 0(e)) or to the curved regions 
(Fig. |2|((i)). So they constitute at this stage the most 
important difference between the present 2D model and 
quasi ID systems of weakly coupled chains. We haven't 
considered Umklapp terms in this work. On physical 
grounds, they are expected to play a very important role, 
in the appearance of Mott insulating phases. The main 
reason for discarding them was to keep the space of effec- 
tive Hamiltonians as small as possible, since even the sim- 
plest situation manifests already a fair amount of com- 
plexity, as we shall illustrate below. Some methodolog- 
ical questions are better addressed in a simple setting, 
which will help in further investigations of more realistic 
models. 

Given the resemblance of our model with quasi ID 
systems, it is natural to recast these effective couplings 
in terms of two species of fermions, which shall be re- 
ferred to as "right" and "left" fermions, as shown on 
Fig. ||(/). This partitioning is to a large extent arbi- 
trary. It is motivated by the desire to consider the two 
parallel Fermi surface segments of a given orientation as 
different species of electrons. Let us consider a situa- 
tion where ki and /c4 are on the "right" vertical seg- 
ment and fc2 and ky, are on the "left" one. According to 
the previous discussion, they are connected by the two 
processes 52(^1, ^2, ^3, fc4)4(^4)c|.,(A:3)cr'(fc2)c^(fci) and 
gi{ki,k2,ki,k:i)c\{k^)cl,{ki)cr'(k2)cr{ki) (see Fig. |). 
In the previous expressions, and in all this article, ex- 
cept otherwise stated, repeated indices are to be summed 
over. c\{k) is of course a fermion creation opera- 
tor, with momentum k and spin projection r. The 
Pauli principle enables us to write the second term as : 
-gi{ki,k2,ki,k:i)cl,{kA)c\.{k'i)cr'{k2)cr{ki). So the 51 
interaction can be cast in a current-current type (right 
goes into right and left goes into left) which reads : 



31(^1,^2,^4,^3) c\{ki)cl,{k:i)Cr'{k2)Cr{ki) 



-E'^AM^A>'cI(fc4)4,(fc3)v(fc2)cM(fci; 



a=l 



(1) 



Here a"^ , a^ and cr" are the usual Pauli matrices. Note 
that the BCS couplings can also be treated in the same 
way. In the limit where all excited states have been in- 
tegrated out, we therefore have the following form of the 



effective Hamiltonian 



^i^t = -^ Yl 4i+fc2,fc3+fc4 X |5'=(fci,A:2,A:3,fc4)c};^,.(fc4)c[_^,(fc3)cL,r'(fc2)cR,r(fcl) (2) 



ki,k2,k3,kA 



In this expression, the wave vectors are chosen in a very 
narrow sheU around the Fermi surface, and the subscripts 
R and L denote the side of the partition, g'^ and g^ will 
be referred to as charge and spin couplings respectively. 
We quantize the system in a finite box of volume fi, so 
the k vectors are chosen on a discrete set, and therefore 
Ski+k2,k3+ki is a Kronecker symbol. With this normaliza- 
tion : {c7-(fc), c^,(fc')} = Sk,k'Sr,T'- In a- renormalization 
group approach, it is common to introduce some cut-offs. 
Since we shall be working at the one-loop approxima- 
tion, we have decided to use a single high energy cut-off 



J 



r 



which will be denoted by A. This means we consider 
only the one-particle states lying at a distance smaller 
than A to the Fermi surface. Following Anderson's "poor 
man scaling" procedure, we will compute the effect of a 
small change in A on the interaction vertex, and show 
this change can be compensated by a small modification 
of the effective couplings. So the effective Hamiltonian 
written above (where states were restricted to the Fermi 
surface) motivates the following form for the finite cut-off 
interacting Hamiltonian : 



(A) 

E 

fel,fc2,'C3,fc4 



^i,!^(A) = ^ E '5fe,+fc„fc3+fe,{gX(^fcl,^fc2,Pfc3,Pfc4)4,,(fc4)cI,,,(fc3)cL,r'(fc2)cR,.(fci) (3) 

+gl(7'fci,7'fc2,7'fc3,7'fc4)a^^a^,^,4 ;,(fc4)c[ ^,(/c3)cL,^'(fc2)cR,^(fci)} . 



In this expression, the wave vectors lie in the strip of 
width 2A around the Fermi surface. The charge and 
spin couplings acquire a A dependence, as physical ob- 
servables (Green's functions) are required not to de- 
pend on A. Following arguments exposed for instance 
in Shankar's review, we shall neglect the dependence of 
g]^ and g\ with respect to the distance of the fc^'s from 
the Fermi surface, since such a dependence is irrelevant 
in the low-energy limit. We have therefore introduced 
the projection V onto the Fermi surface, along normal 
directions. In the same spirit, we shall neglect any fre- 
quency dependence of these functions. Note that the 
BCS couplings are naturally taken into account in this 
description, if we choose ki + k2 = = k^ + k^. 

In a RG treatment, the effective couplings are obtained 
as solutions of a set of differential equations involving 
derivatives with respect to A. The initial conditions of 
these effective couplings' flow have to be specified on the 
basis of a microscopic model. The initial value of A will 
be denoted by Aq, which is then of the order of the inverse 
lattice spacing. As the previous discussion has shown, the 
space of effective couplings is extremely large, since they 
correspond to continuous functions of two or three vari- 



ables. Instead of trying to explore the whole space, we 
shall restrict ourselves to a very simple two-dimensional 
set of initial conditions, namely g\ (fci, fc2, A:3, A:4) = G"^ 
and g^ (fci, A:2, A:3, fc4) — G^. For a local initial interac- 
tion, the microscopic couplings arc indeed independent 
of momenta. The standard Hubbard model corresponds 
to G^ = U/2 and G^ = —U/2, where U is positive (resp. 
negative) in the case of local repulsion (resp. attraction). 
If we were to start with an initial pair potential whose 
Fourier transform will be denoted by U{q), a reasonable 
choice for the initial couplings would be given by : 



gXo(fci,fe,^3,fc4) = U{k4-ki) - \ij{k'i -ki) 
9Ao{ki,k2,k3,k4) = -^U{k3 - ki) 



(4) 



B. ONE-LOOP RENORMALIZATION GROUP 
EQUATIONS 

Because of the very simple form of the effective Hamil- 
tonian, the one-loop corrections to the irreducible two- 
particle vertex functions are given by the sum of the two 
contributions attached to the graphs shown on Fig. 0. 



Examining the particle-particle and particle-hole density 



ki+P 



k\+P 



h+P 



FIG. 3: The two graphs which contribute to the flow of ef- 
fective couplings at the one-loop approximation : (a) corre- 
sponds to the Peierls channel and (fe) to the Cooper chan- 
nel. As often, we shall use full lines for "right" electrons and 
dashed lines for "left" electrons. 



of states, we find that these bubbles are diverging log- 
arithmically as a function of the high energy cut-off A, 
for some special configurations of external leg momenta. 
In the Cooper channel, this occurs if the total incoming 
momentum is vanishing or if it is parallel to one of the 
flat Fermi surface segments. In the particle-hole channel, 
the divergence is obtained if the transferred momentum 
ki — ks is of the form {±2K, ky) or (fca;,±2i^), with \kx\ 
and \ky\ smaller than 2iCtan8. With this type of exter- 
nal configurations, we have the simple result : 
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Ar(fci,u;i, . . . ,^4,0^4) = 
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dfj.{ki +p) 



UJl - W3 



F.s. 27r 
dfi{ki +p) 

F.s. 



gK{ki,k2,k2 -p,ki +p)gA{ki +p,k2 -p,k3,ki) 



2-K 



9K(ki,k^ +p,k3,ki +p)gA{ki + p,k2,k3 + p,k4) . (5) 



r 



In these expressions, the fc^'s are assumed to be on the 
Fermi surface, with fci and ^4 on the "right" side of the 
partition and ^2 and ^3 on the "left" side. In the first 
integral (Cooper channel), p is chosen so that ki +p lies 
on the "right" part of the Fermi surface, and fc2 — p on 
the "left" part. In the second (Peierls channel), ki + p 
again lies on the "right" side, but now k^ + p lies on the 
"left" one. The one-dimensional measure d/i(fci -I- p) is 
just the infinitesimal length element on the Fermi sur- 
face. We have not explicited the charge and spin aspects 



J 



of the couplings at this stage to keep the equations sim- 
ple. Note that if the external momenta are not on a 
single pair of flat segments, ki — k^ is not a nesting vec- 
tor, so only the first term contributes. In fact, for such 
processes, we have ki + k2 = 0, which means they are 
of BCS type, and they renormalize only in the Cooper 
channel. From these expressions, we see the irreducible 
two-particle vertex function remains unchanged (to this 
one- loop approximation), provided the effective couplings 
depend on A according to : 



dgA{ki,k2,k3,k4) 
ainA 



2'kvy 



2lTVp 



d^(fci +p) 



F.s. 
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dnjki +p) 

2tt 



gAiki,h +p,k3,ki +p)gA{ki +p,A;2,^3 +P,k4:) ■ 
I 



(6) 



As previously, all the wave vectors on the RHS which 
appear as arguments in effective couplings are on the 
Fermi surface. It is now important to keep track of the 
charge and spin channels of the interaction. We obtain 
two equations with the same structure as Eq. (O), but 
the quadratic form on the RHS should now be written 
according to Table |. Finally, these equations are still 
difficult to solve, even on a computer. So the integrals 
on the RHS are replaced by finite sums, which amounts to 
discretizing the Fermi surface with a finite set of patches. 
Let us denote by N (resp. M) the number of these 
patches on each straight segment (resp. curved arc) of 
the Fermi surface. We therefore have 2{M+N) patches of 



each type (right and left). If we assume patches of equal 
length, the shape of the Fermi surface, i.e. the angle O, 
is related to M/N by : M/N = (tt - 49)7(4 sin 6). Us- 
ing the "tomographic" parametrization of the couplings 
illustrated on Fig. 0, we may replace the continuous 
functions .gA(fci, ^2, k^, k^) by "vectors" in coupling space 

Cg {I, J). This notation means that a pair of particles 
on patches / on the right and J on the left scatters into 
another pair of patches I -\- S on the right and J + 6 on 
the left. This labelling scheme is consistent with the con- 
straint of momentum conservation. We therefore get the 
discretized version of Eq. (||) : 



TABLE I: Quadratic form on the RHS of Eq. (|) when one 
takes account of the spin. 





First term 
(Cooper Channel) 


Second Term 
(Peierls Channel) 


Flow of 
Charge couphngs 


glgl + 3(7a<7a 


glgl + 35a5a 


Flow of 
Spin couplings 


9lg% + g%9l 
-^glal 


glgl + glgl 



p 



9 In A 167r2vF(M + iV) 



^ /c(/, J, <5, «)Ci^)(/, J)&t\{I + a,J + a) 



J2 /P (I, J, S, a)Ci^^ (/, J + <5 - a)C^ti (I + «> •^) I 

a ) 



(7) 



Here, P is the perimeter of the Fermi surface, so P = 
[(27r — 80)/ cos 9 + 8 tan Q]K. /c and /p enforce phase 
space restrictions for the Cooper and Peierls channels. 
In the Peierls channel, our previous discussion has shown 
that patches I, I + a and 1 + 6 have to lie on a given flat 
segment on the right side, with patches J, J + S — a and 
J+^ on the opposite parallel segment. If this condition is 
met, /p is set equal to one, and vanishes otherwise. In the 
Cooper channel, /c is equal to one if a similar condition 
is satisfied. But this channel allows for processes which 
connect one set of parallel flat segments to the remainder 
of the Fermi surface. In this case, I = J and /c is also 
equal to one. In any other situation, /p vanishes. By 
using Table | and with a minor change in the couplings' 
normalization, we obtain the flow equations (2) and (3) 
of our previous workll3, which we won't reproduce here 
to save space. 

We have performed numerical evaluations of the dis- 
cretized equations for various initial conditions of the 
form : 




Cf^°\l,j) 



C^ and C: 



i(Ao) 



(/, J) = C' 



(8) 



FIG. 4: A "tomographic" labelling of Fermi surface patches 
in the A'^ = 3 and A/ = 1 case. In this description, two op- 
posite points on the Fermi surface are assigned the same pos- 
itive integer, but they are distinguished by their "right" (R) 
or "left"(L) nature. Arrows illustrate the Cs{I,J) = Ci(l,2) 
coupling (see text). 



Again the simple Hubbard model corresponds to C"^ + 
C^ = Q with C^ positive in the case of repulsion. In 
most cases we find that the effective couplings diverge 
after reaching a low-energy scale Ac < Aq. Although 
we are working in a zero temperature formalism, it is 
natural (and often done in the literature) to interpret 
Ac as a typical temperature scale for the onset of some 
kind of long range order, either in the particle-particle 
channel (superconductivity) or in the particle-hole chan- 
nel (charge or spin density waves). As will be described 
in detail in Sec. VB below, the response functions for 



the corresponding order parameters can be evaluated in 
the one-loop approximation, and are found to diverge in 
general as power laws of ln(A/Ac), where A is a typical 
frequency scale at which the system is probed. 



In order to get a better understanding of these one-loop 
RG flows, we have found useful to decompose them in two 
different (although not unrelated) aspects : the evolution 
of the magnitude of the couplings, and the evolution of 
the direction in the space of all possible couphngs. This 
is possible because one-loop flows have a nice invariance 
under dilations of the effective couplings. To see this, 
we note that our RG equations take the simple fornr : 
dtgi = Aijkgjgk- Here, the flow's time t is related to 
the current cut-off A by t = ln(Ao/A). The giS are for 
instance the set of (7^(7, J) 's defined above and used in 
numerical computations. The initial conditions (micro- 
scopic couplings) correspond to i = 0, and t increases as 
the high energy cut-off is gradually reduced. Since the 
RHS of the one-loop RG equations is homogeneous and 



quadratic in the couplings, we have the simple property 
that if gi{t) is a solution so is gl (t) = Xgi{Xt) for any A. 

This implies that g| {t)/g, (t) = gi{Xt)/gj{Xt), namely 
that the direction of the coupling vector follows a trajec- 
tory which is independent of the magnitude of the initial 
couplings. It is therefore natural to introduce normalized 
couplings hi{t) defined by 



Kit) 



N{tr 



withAA(t)= \Y,9^{tf 



(9) 



We then have for any positive A, and with obvious nota- 
tion : 7V(^)(0 = XN{Xt) and hf\t) = hi{Xt). Next we 
wish to write down differential equations for A/'(i) and 
for hi{t) instead of gi{t). It is straightforward to show 
that one has : 



dghi — Aijkhjhk — {Ajkihjhkhi)hi 
dsAf ^ {A.jkihjhkhi)M 



(10) 



with s being the "time" adapted to the flow of hi, de- 
fined by As = Af{t)dt and s{t — 0) == 0. Changing 
from i to s amounts to zooming near the singularity at 
ic = ln(Ao/Ac), and s can therefore be taken as large as 
possible without actually reaching the divergence of J\f. 
However, we should keep in mind that our RG equations 
are valid only if Af{t) is not too large. How does this 
translate on the flows regarded as functions of s ? From 
Eq. (p^), we have : 

Af{s)=Af{0)e^^'\ with 

F{s) = f ds'Ajkihj{s')hk{s')hi{s') . (11) 

In most cases, F(s) becomes infinite as s goes to infinity, 
so there will be a typical value Smax of s beyond which 
the one-loop equations are no longer valid. But Smax 
depends on the magnitude of the initial couplings, and 
can be made as large as desired, if Af{Q) is replaced by 
A7V(0) with A small enough, while keeping the initial di- 
rection hi{0) fixed. Part of this work deals with the large 
s limit, which may be equivalently viewed as a very weak 
coupling limit. 

Since the flow of the normalized couplings defined 
above is constrained to stay on a compact space (the unit 
sphere), it is expected to reach either limit cycles or fixed 
points in the large s limit. We have never observed cycles, 
but fixed points (in the space of directions) are always 
found for large s. These correspond to solutions of the 
type g'i,{t) = u.J{tc-t), with u^ satisfying Ui = AijkUjUk- 
Fixed directions are interesting since they are character- 
ized by a well defined set of exponents for the growth 
of various response functions in both Cooper and Peierls 
channels (see Sec. VB ). This pointi-aLview has already 
been emphasized by several authorsJlJO'E^I The new fea- 
ture of the present work lies in our attempt to reach 
conclusions for large values of the number N of patches. 



as required for the description of a macroscopic 2D sys- 
tem. Note that when the system is flowing along a 
fixed direction, F(s) defined above is equal to s/c, with 

c^{J2^u^'^t^^■ SoAA(s) =AA(0)exp(s/c). We also have 
t{s) = c[l - exp(-s/c)]/7V(0), so t^ = c/Af{0). 



C. NUMERICAL RESULTS FOR THE 2D 
MODEL WITH NESTED FERMI SURFACE 

We have performed numerical evaluations of the RG 
flows for a large set of ratios C^/C^ of the initial cou- 
plings and for various values of M and N, where M /N 
is the length of the curved regions divided by the length 
of the flat regions on the Fermi surface. We have used 
the set of Eqs. |lO| for the normalized coupHngs hi. Be- 
sides the theoretical motivations just discussed for intro- 
ducing them, they have the practical advantage of re- 
maining bounded during the RG flow, by contrast to the 
giS which diverge at the finite "time" t^- In Fig. ^, we 
show the corresponding phase diagram for TV = 6 and 
M = 0, together with typical flow patterns (using the s 
variable and normalized couplings) for the charge density 
wave (Fig. ||) and for the d-wave superconducting phase 
(Fig. m). These results have already been discussed. t£l 
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FIG. 5: Phase diagram for a nested system . The flows of the 
couplings for the charge density wave phase and the d-wave 
superconducting phase are shown in Figs. H and M respectively. 

In this work we shall focus on the two regions which oc- 
cur for positive charge couplings, and when the absolute 
value of C^/C^ is not too large. In these conditions, the 
BCS type couplings (i.e. with a vanishing total incom- 
ing momentum) are strongly reduced in the early stages 
of the flow. The couplings which are building up are 
of the Peierls-type, meaning they can be interpreted as 
processes with a transferred momentum Q — {±2K, 0) 
or Q — (0,±2A') connecting one flat region to another. 
Note that the final behavior of the flow (for very large s) 
looks very different according to whether the initial spin 
coupling C^ is positive or not. If C^ > 0, and C/C"^ not 
too large, the system flows towards a charge density wave 
instability. As will be discussed below, the approach to 
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FIG. 6: Flow of all the charge (top) and spin (bottom) cou- 
plings, in the case of the charge density wave phase of Fig. p| 
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FIG. 7; Flow of all the charge (top) and spin (bottom) cou- 
plings, in the case of the d-wave superconducting phase dia- 
gram of Fig. 0. 



the transition is qualitatively well described by a RPA 
like description which singles out the Peierls channel. If 
C < 0, and IC^/C^l not too large, we observe three 
stages in the RG flow, instead of two for C^ > 0. Af- 
ter a transient regime which eliminates most of the BCS 
couplings, we find a rather long intermediate regime in 
which the dominant correlations are spin density wave 
like. This regimeJs very similar to the one obtained by 
Zheleznyak et al.tj In agreement with them, we find a 
tendency for the dominant couplings to be larger for in- 
coming and outgoing electrons close to the end points 
of the flat regions. But running the flow closer to the 
singularity shows that this regime is eventually unstable 
towards the formation of a very different type of correla- 
tions. Analysis of the corresponding couplings shows that 
events involving all the incoming and outgoing particles 
on the end points of the flat regions dominate all others. 
Some of these processes are illustrated on Fig. g. So the 
stable fixed direction reached by the flow in this regime 
corresponds to an eight-patch model, independently of 
the actual number of patches involved in the numerical 
computation of the flow ! In physical terms, these eight 
patches are reminiscent of the "hot spots" observed in 
angle-resolved photoemission spectra of high tempera- 
ture superconductors, for electron wave vectors in the 
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FIG. 8: Some of the few processes which dominate in the final 
stage of the flow (d-wave superconductivity fixed direction) 
for C"= < 0, C" > and |C7C"| not too large. 



vicinity of the Eermi surface intersections with Brillouin 
zone boudaries.E3 This final regime is also remarkable 
since the leading instability is in the Cooper channel, 
with the formation of singlet particle pairs of d-wave like 
orbital symmetry. So we interpreted the flow in this part 
of the phase diagram (in which the usual Hubbard model 
with local repulsion falls) as manifesting a cross-over phe- 
nomenon : between a metallic high temperature regime 
and a d-wave superconducting phase lies an intermedi- 
ate state with dominant itinerant antiferromagnetic cor- 
relations. Note that superconductivity occurs here in a 
perfectly nest£d|^sjbem. This is to be contrasted with 
other studies, Baaij'B in which superconductivity is pos- 
sible only when the nesting is not perfect. With an im- 
perfect nesting, contributions from the Peierls channel 
are suppressed at low energies, where only the Cooper 
channel survives. 

One of the motivations for this work was the desire to 
understand to which extent such a picture has a chance 
to apply to real macroscopic 2D systems. Three types 
of objections have been considered. The first one is in- 
spired by the analogy just mentionned with some as- 
pects of high temperature superconductors. One of the 
most striking features of this high-Tc cuprates is the pres- 
ence of the strange pseudo-gap phase in the underdoped 
compounds. At least two types of scenarios have been 
proposed to describe this regime. In the most popular 
one, local Cooper pairs are formed around a temper- 



ature T* above Tc, as suggested by a clear reduction 
of angle resolved photoemission intensity at the chem- 
ical potential energy, near the Brillouin zone boundary, 
where the single-particle spectral weight is pushed below 
the Fermi level. Tc is then interpreted as^he onset of 
phase coherence of these preformed pairs .EZl In another 
scenaricS superconductivity is argued to compete with 
another instability, whose typical onset energy scale van- 
ishes as one goes from the underdoped to the overdoped 
regime. Our intermediate regime would fall in this latter 
type of scenario. However the simple model considered 
here doesn't seem to generate a spin gap in the inter- 
mediate regime. Our numerical results indicate that the 
static spin susceptibility at the nesting vector Q steadily 
increases as the typical energy scale is reduced, showing 
no evidence of a spin gap. The cross-over to a d-wave su- 
perconductor at lower energy is simply deduced from the 
fact that the superconducting response function grows 
faster than the Neel one as s goes to infinity. Clearly, 
since we are dealing with a low-dimensional system, we 
have to investigate the effect of long wavelength fluctua- 
tions of the Neel o rder parameter. This question is briefly 
addressed in Sec. VC below, and we find, within a sim- 



ple version of the intermediate fixed direction exhibiting 
SU(Af) symmetry, that this "Goldstone mode" fluctua- 
tion induced gap vanishes like exp(— iV) as N goes to 
infinity. This is an interesting (though not very surpris- 
ing) conclusion : other physical ingredients are required 
to explain the pseudo-gap regime of high-Tc cuprates. 
A natural candidate for this are Umklapp processes 
discussed within a RG approach by Honerkamp et al 
On intuitive grounds, the spin gap corresponds to short 
range singlet correlations. Since there is no frustration 
on a square lattice, the strongest source of fluctuations 
comes from the motion of holes which do a lot of damage 
in an ordered magnetic state. Clearly, this picture makes 
sense for a moderately doped Mott insulator, and Umk- 
lapps are the natural way to describe the vicinity of the 
Mott insulator in a weak coupling RG approach. Umk- 
lapps have already been considered by Zheleznyak et al, 
for a very similar situation, and were shown not to affect 
much the spin density wave like regime, in the case of a 
half-filled Hubbard model. These authors haven't con- 
sidered doped situations where the Fermi surface would 
remain perfectly nested, but with twice the nesting vec- 
tors Q not exactly on the reciprocal lattice. It would be 
very interesting to study whether a pseudo-gap behavior 
appears in spin or uniform charge responses for such a 
case. We leave this for future investigation. 

Leaving aside the question of whether we may obtain 
an intermediate regime with a pseudo-gap, we return 
to the question of assessing the possibility of observing 
cross-over phenomena involving an intermediate direc- 
tion as described above. As will be shown below, there 
does not seem to be anyj^roblem with the numerics on 
which our previous workllj is based. This would be our 
second objection. Indeed, the example of an intermedi- 
ate direction with SU(iV) symmetry is treated in great 



detail below. In the large N limit, we have established 
analytically that this system admits {N— 1) relevant per- 
turbations for the flow of the normalized couplings. So 
the existence of these instabilities and the tendency for 
the system to keep only a small finite set of "hot spots" in 
the infinite s limit seem firmly established for any finite 
value of N. 

The most serious objection, which hasn't been ad- 
dressed before, deals with the infinite N limit. Recall 
that patches were introduced to replace integrals of func- 
tions of momenta chosen on a continuous Fermi surface 
by finite sums. The validity of such an approximation 
depends on the way effective couplings vary with exter- 
nal momenta. But the large s regime exhibits a very 
singular variation, since a hot spot picture emerges. In 
this situation, it is highly desirable to keep N as large 
as possible. As will be shown, the relevant perturba- 
tions around SU(A^) fixed direction of spin density wave 
type become marginal in the infinite N limit. This has 
the consequence that the RG time required to leave the 
vicinity of the intermediate fixed direction diverges as 
N goes to infinity. According to the arguments detailed 
above, this implies we need to restrict ourselves to van- 
ishingly small initial couplings as N becomes infinite, in 
order to observe the final stage of the fiow in the validity 
range of one-loop RG equations. At the present stage, 
we don't have any strong argument to counteract this 
objection. So the present work reinforces the validity of 
the more traditional interpretation already presented by 
Zheleznyak et al : in a real system, the instability of the 
couplings occurs while they are close to an intermedi- 
ate fixed direction. The remaining parts of this article 
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FIG. 9: The dominant couplings A{I), C{I,J) and D{I,J) 
involving the vertical Fermi segments in the regimes of density 
wave like instabilities. The patch indices / and J refer now 
to the projection of the momentum along the vertical axis. 
Two Fermi surface points linked by a dashed line thus have 
the same patch index. 

are mostly dedicated to the physical properties of these 
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intermediate fixed directions. 

Let us now motivate the strategy we have foUowed. 
A detailed analysis of the effective couplings in the inter- 
mediate regime shows that the BCS couplings which con- 
nect a pair of parallel segments of the Fermi surface to 
its complementary parts are extremely small compared 
to the dominant couplings. For a system with N = 7 
and M = 0, we found a typical ratio of 10~^ between 
these two classes of couplings. So, we shall adopt the 
approximation already made by Zheleznyak et al, i.e. we 
shall focus our attention on a given pair of paradlel seg- 
ments, thus restricting ourselves to a ID setting .Il3 Next, 
we find that the dominant couplings in this regime are of 
the three types depicted on Fig. g. The D{I, J) couplings 
are reinterpreted (using the Pauli principle) as scattering 
events with a momentum transfer Q — (±27^,0). Note 
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FIG. 10: Example of a typical ID RG flow pattern for A'' = 
16, and initial condition G" — —C. Only charge couplings 
are represented here (see Fig. h]l for spin couplings). The 
four figures show the behavior of D, C, A, and non-principal 
couplings respectively. We have represented the couplings 
that survive in the end of the flow (that is involving patch 
numbers 1 and 16) with triangles. With circles and squares, 
we have shown the behavior of two charge couplings that are 
a sort of "envelope" of the couplings. By this we mean that 
lots of couplings are found between the values of these two 
couplings, as was also the case in Fig s, o and M where all 
couplings were shown (see also Fig. n4 for a detailed view of 
the distribution of the values of the couplings). 
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FIG. 11: Same as Fig. ffd for the spin couplings. 



that a symmetrical pattern is obtained for the horizontal 
segments. For the sake of simplicity, it is not shown on 
the figure. This classification of the couplings according 
to their magnitude in the intermediate regime of the flow 
is well illustrated on Figs. HG to O. These curves have 
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FIG. 12: Zoom of Fig. |l| for times s between and 200. 
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FIG. 13: Zoom of Fig. for times s between and 200. 

been obtained in a purely ID setting, since we haven't 
found much difference between the full 2D case and this 
simplified version which involves N coupled ID chains. 
The corresponding formalism is outlined in section III 
below. Roughly speaking, the initial transient of the flow 
is observed for s < 50, and the intermediate regime for 
50 < 1000. The selection of A, C and D couphngs with 
respect to the remaining ones (which will be referred to 
as non-principal couplings) occurs clearly for s between 
50 and 100 (see Figs. |l| and |l|). 

To understand better the pattern of the A{I)'s, 
C{I, J)'s and D{I, J)'s couphngs, we notice that for most 
values of / and J, the sign of these couplings, in both 
charge and spin sector are well defined, and that the mag- 
nitudes do not vary strongly in a given sector, as shown 
on Figs. O and U3. So in the sequel wc shall consider 
fixed directions with A, C and D couplings chosen not to 
depend on / and J . As will be shown in Sec. V A below, 
this leads us to situations where A'^'^ — C^'^ + D'^'^. If 
such a relation holds, the effective Hamiltonian exhibits 
a large SU(iV) symmetry in the space of "patches" . This 
phenomenon of enlarged symmetries is rather common 
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in one-dimensional quantum field theories. A striking 
example is the appearance of a SO (8) symmetry in a two 
channel model which may be relevant to the descriptioa. 
of low-energy properties of metallic carbon nanotubes.t£l 
Most of the sequel will be dedicated to these special 
points for which a large amount of information may be 
obtained by analytical calculations, at least in the large 
N limit. In Fig. |lj, we compare the actual distribu- 
tion of dominant couplings of C and D type (both for 
charge and spin sectors) in the intermediate regime (for 
s — 100) with the expected values under the assump- 
tion of SU(Af) symmetry. Clearly the actual flow doesn't 
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FIG. 14: Normalized repartition of D and C, charge and spin 
couplings in Fig. [13 and Fig. ISlat time s — 100. Each bar rep- 
resents the normalized number of couplings taking the value 
theywould have on the SU(A'') (-1-,-) fixed direction (see Ta- 
ble O). The fine vertical line represent the mean A^, and the 
arrows represent the variance. 

reach exactly an SU(Af) symmetrical fixed direction, but 
the histograms of C and D couplings are well centered 
around the predicted value for SU(iV) with N = 16. As 
we shall see, the Luttinger type couplings C^'^ are of or- 
der 1/A^ whereas A'^'^ and V^'^ are of order unity if TV 
is large. This remark suggests another way to analyze 
the intermediate fixed direction. If the C(/, J) are as- 
sumed to vanish, the D(/, J) couplings are renormalized 
in the Peierls channel only. So an RPA type approxima- 
tion should give a rather good picture for the vicinity of 
these fixed directions. To test this idea, we present in 
Fig. ^ a comparison between the average contributions 
of the Cooper and the Peierls channels to the flow equa- 
tions, as a function of the new "time" variable s. More 
precisely, we define for each channel (charge c or spin 
s) a ratio (3 measuring the relative importance between 
Cooper and Peierls channels as follows. We first rewrite 
the flow equations in the charge (spin) channel, distin- 
guishing between the Peierls and Cooper contributions, 

in the form : digf""^ = ^yfe^'^5j5fe - A^^^l^^'^jgu- With 
obvious notations, this is cast in a vectorial form as : 
Qgcis) ^ ^c,c(s)(-^^ _ ^P,c(s)(g-)^ Finally in each chan- 
nel we define /?'=(") = || A^''=(")||/||AP^'=(')||, where ||.|| is 
the Euclidean norm. Clearly, the Peierls channel domi- 
nates the flow as soon as s > 50 and for s < 1200 (see 
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FIG. 15: Evolution of fi"^ and (3'^ (see text) with time, for 
different initial conditions, and for N = 16. (a) and (b) : 
C^ = — C = 0.1 (6) is a zoom on small times of (a). At 
s — 0, P" = oo. This is peculiar to the initial condition for 
which C^ — —G^ . In (c) we show what happens for C^ = 
—G" /2 = 0.1 at s = 0. Finally in (d) we show the flow for 
initial condition G"" = C = 0.1 



Fig. |l^(6)). One can numerically check that in the in- 
termediate regime, the values of the /? ratios scale as 
1/iV when one varies N. In this situation, neglecting 
the Cooper terms yields a smooth manifold of possible 
flxed directions in the space of D{I, J) couplings, as will 
be shown below in Sec. IV. This manifold contains the 



special SU(A'^) symmetry points discussed above as spe- 
cial isolated points. So, the RPA picture shows that 
the SU(Af) symmetry doesn't really emerge (see again 
Fig. |l^). But these special points are nevertheless very 
interesting, since the actual flows get close to them in the 
intermediate regime, and their underlying symmetry is a 
source of tremendous analytical simpliflcation. 

Before going further, we would like to stress the 
following importaat point, actually already made by 
Zheleznyak et al.Ej Although the intermediate regime is 
very well described by a RPA picture, it is quantitatively 
wrong to neglect the Cooper channel in the early stages 
of the flow, as indicated by Fig. |l5|. This would give 
a higher value of Ac than obtained from the complete 
"parquet" solution of the one-loop equations. But also, 
as seen when comparing Figs. ^ and |l^ with Figs. |l^ 
and n3, the initial transient has to be accurately treated 
in order to determine the actual fixed direction chosen 
by the system in its RPA like intermediate regime. The 
interference between the Cooper and Peierls channels is 
indeed responsible for the observed splitting among the 
values of the dominant charge and spin couplings in this 
regime. A simple RPA (i.e. no Cooper terms at all) start- 
ing from fc-independent couplings preserves this constant 
character of A and D couplings along the flow, by con- 
trast to what happens in the "real" one-loop flow. 

To summarize, we shall follow two complementary ap- 
proaches. In Sec. IV we shall flrst set up a perturba- 
tion analysis around the RPA in the Peierls channel, the 
Cooper channel being considered as a small perturbation. 
In this picture, the RPA exhibits a continuous manifold 



12 



0.1 

^^0.05 



0.1 

^" 0.05 





1 

a 

y 

1 


1 


1 


1 


1 


50 


100 

s 


150 20 


1 

n 

y 

1 


1 


1 


1 


1 





0.04 
GSpO.G2 




50 100 150 200 







50 100 150 200 




FIG. 16: Same as Fig. [L3in the RPA. We have not represented 
the flow up to s = 1600 as in Fig. hd because the couplings 
do not evolve any more. It is remarkable that all D and A 
couplings remain equal. 
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FIG. 17: Same as Fig. h^ in the RPA. 



of fixed directions, which may be destabilized by the weak 
Cooper terms, but this final stage involves very long RG 
time scales. Alternatively, we may say that the system in 
the intermediate regime reaches a neighborhood of a fixed 
direction (for the full RG flow this time) with SU(7V) 
symmetry. This direction admits a set of {N — 1) rele- 
vant perturbations which become marginal in the infinite 
N limit. This will be the object of Sees. |^ and 0. We 
shall first of all set up our notations in the ID case and 
study in detail its set of fixed points in the usual sense. 



III. ONE-LOOP RG FOR N COUPLED CHAINS 

A. GENERAL SETTING AND RG EQUATIONS 

We have seen that pairs of nested flat regions on 
the Fermi surface produce a complex dependence of 
the electron interaction as a function of the typical en- 
ergy scale. These phenomena arise from the mutual 
couplings between the particle-particle and particle-hole 
channels. These main features are already present in one- 
dimensional systems, for which RG equations are a little 
simpler. Furthermore, the previous discussion of our nu- 
merical results suggests that many couplings become neg- 



ligible in comparison to a much smaller set of dominant 
couplings, as the system flows towards lower energies. 
We have seen that often these dominant couplings corre- 
spond at least approximately to a very large SU(A'^) sym- 
metry of the effective interaction, at some intermediate 
stage along the flow. Keeping only these couplings brings 
dramatic simpliflcation, especially in the one-dimensional 
case of N coupled chains, where analytic approaches be- 
come available. 

To be specific, we shall concentrate on systems of N 
chains without transverse hopping (to generate a flat 
Fermi surface) and coupled by two-particle interactions 
between right and left movers, with a small transferred 
momentum along the chain directions. Each chain may 
be seen as a device to describe the electronic states in 
the vicinity of a point chosen on a fiat section of the 
Fermi surface. In what follows, the chain direction is un- 
derstood to be perpendicular to the Fermi surface. The 
kinetic energy is taken to be iJo, with : 



-ffo = WF ^ |(fc - fcF)cR,j,^(fc)cRJ,<T(fc) 

-(fc + fcF)c^,(fc)cL,/..(fc)}. (12) 



k,I,a 



Here up and fcp stand for the Fermi velocity and wave- 
vector. The operators c}g^n)ia-^k) and cr{,(L),/,cr(fc) are 
usual creation and destruction operators for right (left) 
moving spin 1/2 fermions. k is the momentum along 
the chains, / is a "color" index which refers to a given 
patch on the Fermi surface (according to the previous 
section) and a is the z-component of the spin. The sum 
over k runs for k = 21:7x1/ L, m being any integer, and L 
denoting the length of the chains. We interpret / as the 
momentum in the transverse direction of the chains. It 
is chosen in the set {1, . . . , N} (see Fig. |l^). The most 
general spin rotation invariant two-particle interaction 
Hamiltonian, involving scattering events between right 
and left movers may be written as : 



-ffint — 
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+GUl,J)Sn,iAQ)SL,jM-<l)} ■ (13) 

In the summation, both / and J belong to {1, . . . , N} 
and S is chosen in such a way that both I+S and J—S be- 
long to this set of patches (see the illustration on Fig. ^ . 
Note that this way of writing the two-particle interaction 
enforces momentum conservation, in both perpendicular 
(q variable) and parallel {S variable) directions with re- 
spect to the Fermi surface. The restrictions on I, J, I + S 
and J — S imply that we are ruling out Umklapp pro- 
cesses along the Fermi surface. Therefore this model is 
not relevant for quasi one-dimensional conductors which 
also exhibit disconnected Fermi surfaces with nearly fiat 
curves running all across the first Brillouin zone. Again, 
our chains are introduced to describe processes in the 
vicinity of the flat regions of a connected Fermi surface. 
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FIG. 18: Labelling of the patches in the ID case for N — 8, 
according to the interpretation in terms of chains. Arrows 
illustrate the Gs{I, J) ~ G2(3, 8) couphng (see text). 
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The functions G^{I,J) and Gl{I,J) will be referred to 
as charge and spin couplings. With our normalization of 
ifint, these couplings are dimensionless. In this Hamil- 
tonian, we have neglected processes which involve all in- 
coming and outgoing particles on the same branch since 
they don't appear in the flow of the couplings at the 
one-loop level. They commute with Hq, so their main ef- 
fect is to lift some degeneracies between various excited 
states of Hq which often amounts to changing Fermi ve- 
locities. The renormalization of these velocities cannot 
be addressed within a one-loop approximation. So for 



i + s 



I + S 



This is the reason we introduce this kind of open bound- 
ary conditions at the end points of the flat regions, which 
correspond to chains labelled by 1 or N. The quantities 
PRj.siq) and Sjij^siq) are generalized charge and spin 
right currents (and similarly for the L side) defined by : 

fe,r 
SllAl) = E 4„/+5,.(fc + '?)<r'CR,7,.'(fc),(15) 

where a = x,y, z and cr° are the usual Pauli matrices. 
I 



FIG. 19: The particle-hole and particle-particle graphs that 
contribute to the ID RG equations. We have not indicated 
any spin indices. 



simplicity, we assume all our Fermi velocities to be equal 
and constant along the RG flow. We are also discard- 
ing Umklapp processes along the chains, which play a 
dramatic role only for a discrete set of band fillings. 

The total Hamiltonian H — Hq + Hint describes a 
renormalizable field theory. The one-loop flow equations 
are obtained by simply computing the two graphs shown 
on Fig. M, and we get : 



' dtC^il, J) = ^ Ea{Gl{I. J + a- S)G's-M + a,J) + 3G^(/, J + a- <5)G|_J/ + a, J) 

-GUI, J)G-s-aiI + a,J~a)- 3G«„(/, J)Gl^^iI + a, J - a)} 

dtGlil, J) = ^ Ea{'2Gl{I, J + a- S)Gl_^{I + a,J) + Gl{I,J + a- S)G's_JI + a, J) 

+Gl{I, J + a- 5)Gl_JI + a,J)+ 2G^(/, J)G|_„(/ + a,J-a) 

-GUI, J)G'i_Ul + a,J-a)- Gl{I, J)Gl^JI + a, J - a)} 



(16) 



We have set G = 2ttvfG in order to get rid of factors 
2ttvf. As before the dt symbol means the derivative with 
respect to i = ln(Ao/i^), where Aq is a "bare" (i.e. mi- 
croscopic) cut-off, and i^ < Aq is the typical energy scale 
of the process of interaction. In the sequel, we shall also 
analyze the simpler case of spinless fermions, which is 
simply obtained by taking G|(/, J) = in Eq. ([l6|). In 
this case, G^{I, J) will simply be written as Gs{I, J), for 



which the one-loop fiow reads : 

dtGsil, J) = ^ Y.{gM, J + a- 6)Gs-M + «, J) 

-G„(/, J)Gs-a{I + a, J - a)} . (17) 

Note that the dominant couplings in the intermediate 
regime introduced in Sec. |II C| correspond in this language 
to : A{I) = Go(/, J), G(/7J) = Go(/, J) and D{I, J) = 
Gj^i{I,J). It is easy to check they form a closed set 
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under Eq. ( |T^ ) 



B. FIXED POINT SOLUTIONS AND THEIR 
STABILITY 

1. "k-space" Luttinger liquids 

As usual, it is convenient to begin the study of these 
flows by searching for fixed points (i.e. true fixed points, 
not fixed directions). This discussion is logically inde- 
pendent of the main line of this work, but we include it 
here, since it is a rather simple and natural question to 
consider. Furthermore it illustrates, in a different man- 
ner, some difficulties which arise in taking the large N 
limit (see below). In general, we have found a continuous 
family of fixed points which may be seen as generalized 
Luttinger liquids. They are characterized by the property 
that all couplings for which 6^0 vanish. We remind that 
6 denotes the transferred momentum along the Fermi sur- 
face. Furthermore, the spin couplings Gq*{I, J) for 5 = 
all vanish (the star characterizes the fixed points). So 
these fixed points are parametrized by a two-parameter 
family of couplings Gq*(/, J) (still for 6 = 0). Physi- 
cally, they are Luttinger liquids in "fc-space" , in close re- 
semblance with P.W. Anderson's idea of what he called 
a "tomographic" Luttinger liquidJ23 The corresponding 
fixed point Hamiltonian may be cast as : 



H = 



27rwF 



Yl {pRj(5')PRj(-'?) + PL,/(-g)pLj(g)| 
-^ J2 G^o*iI,J)PnAl)PL,ji-q)- (18) 



Here, we have used the notation pvij{q) for pr,/, 5=0(9) 
and similarly for left movers. Since it is quadratic in den- 
sity modes, it is easily diagonalized, yielding a decoupled 
set of N branches of non-interacting sound-like bosonic 
modes. Usual bosonization techniques may be used to 
obtain exact low-energy asymptotic expressions for the 
corresponding fermion propagator and various correla- 
tion functions. But we shall not develop this here. 

Let us now investigate the linear stability of these gen- 
eralized Luttinger liquid fixed points. To do this, we 
parametrize the set of couplings as : 



G3(/,J)-GS*(/,J)+5S(/,J) 
G^,(/,J) = <?,^(/,J), iorS^O 
G|(/,J) = <?^(/,J), for any <5 



(19) 



Here g'^'^(I,J) are assumed to be a small perturbation. 
Taking the flow equations and keeping only the linear 
terms in g variables gives : 



dtg^oiI,J)^Oig^) 

dtgl{I,J) = 0{g^) 

dtg's{I,J) = -wB*{I,J,S)gliI,J) 

dtgliI,J) = -^B*{I,J,5)gliI,J) 



(20) 



where : 

B*{I,J,S)=GriI,J)-Gr{I,J-6) 

-Gg*(/ + <5,J)+GS*(/ + 5,J-5). (21) 

So in order to be stable with respect to small perturba- 
tions, these fixed points are required to satisfy the fol- 
lowing inequalities : 



B*{I,J,6) >0 for any (5 



(22) 



It is easy to show that this condition is satisfied for 
any 6 provided it is satisfied for 6 = 1. Note that the 
perturbations at zero momentum transfer gg'^^I, J) re- 
main marginal around these Luttinger fixed points. If 
we examine their flow in more detail, we find on the one 
hand that 9t (?§(/, J) involves only irrelevant couplings 
ga{P, Q) with a 7^ 0. Therefore dtgQ{I, J) is expected to 
be small, provided the bare couplings satisfy the stability 
condition stated above. On the other hand, we obtain : 
dtgQ{I,J) = {4:/N)gQ{I,J) + irrelevant terms. Assum- 
ing the irrelevant terms vanish as the energy cut-off goes 
to zero, we see that the flow remains in the attraction do- 
main of the generalized Luttinger liquid provided all the 
spin couplings remain negative. Otherwise, some of these 
couplings are expected to become marginally relevant, 
with the opening of a spin gap as a likely consequence. 

The main result of this discussion is the possibility (for 
any finite N) to stabilize generalized Luttinger liquids 
in "fc-space" , characterized by processes involving mo- 
mentum transfers exclusively perpendicular to the Fermi 
surface. These fixed points are parametrized by the in- 
teraction function Gq*{I,J) which plays the role of the 
Landau parameters for the Fermi liquid fixed point. Sta- 
bility is ensured by the not too stringent Eq. (]2^), re- 
stricted to the case 5=1. The main difficulty is that our 
results are valid for a finite number TV of Fermi points. 
For a realistic system, we would like to take the N ^ 00 
limit. But this is tricky, as indicated by the 1/N factor 
in the RHS of the linearized RG equations (20). In this 
limit, the ^ 7^ couplings become marginal instead of 
being irrelevant as in the finite N case. In the infinite 
N hmit, Gs*{I,J) is replaced by a function Gd*{u,v) 
of three continuous variables, where u and v denote the 
transverse momenta of the incoming particles, and d the 
transverse component (i.e. parallel to the fiat Fermi sur- 
face) of the momentum transfer. The RG equations take 
the form : 

5tGJi(u,v) = / da|GJj(w, u-fa-(i)G^„„(u-|-a,u) 

-Gl{u,v)G'i,_^{u + a,v^a)]. (23) 

To simplify this discussion we assume the spin couplings 
are vanishing. The tomographic Luttinger liquids dis- 
cussed above correspond to GJ^*(m,i') = 6[d)G*{u,v), 
and the stability condition becomes : 



92 G* 
dudv 



{u,v) < 



(24) 



But clearly these are singular couplings. The main ques- 
tion which has to be addressed is the following : suppose 
the bare (microscopic) couplings are smooth functions of 
the three variables d, u and w ; is it possible to build up a 
5{d) singularity in a finite RG time ? Although wc can't 
make any rigourous statement, our experience with nu- 
merical solutions of RG flows, attempting to increase N 
as far as A^ = 20, indicates that these singularities do not 
appear, unless they are already present in the bare cou- 
plings. So we are tempted to believe that tomographic 
Luttinger liquids do not appear as fixed points for sys- 
tems with smooth microscopic interactions in fc-space, 
that is with short range interactions in real space. This 
example shows clearly that a lot of caution is required 
in extrapolating numerical results obtained for a finite 
iV to the continuous Fermi surface relevant to real life 
situations. 



2. Comparison with "real space" Luttinger liquids 

We would like for completeness to make contact with 
the physics of strongly anisotropic conductors. In the 
previous discussion, chains are seen as a mathematical 
device to discretize a continuous Fermi surface of a two- 
dimensional conductor along its fiat regions. Let us now 
consider a one-dimensional system made of a series of 
N parallel identical chains, and let us assume periodic 
boundary conditions in the perpendicular direction with 



J 
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respect to the chains. To simplify, we assume there is no 
single-particle hopping from one chain to the next. As a 
result, we obtain two flat Fermi surface segments along 
the y-axis (if x is the chain direction) running along the 
whole Brillouin zone, from ky — —n to ky — tt (we choose 
the unit length to be the separation between two neigh- 
boring chains). We consider the following interaction : 



Hiat = / dx > F{m',m;n',n) 
JO , I 

X : V'R,™'(a;)V'R,m(a;)V'L,n'(^)^L,n(a;) : , (25) 

where : O : denotes normal ordering of opertator O. 
ipR,mix) is a fermion destruction operator at point x, 
on chain labelled by the integer m. To simplify this dis- 
cussion, we shall omit the spin degree of freedom. The 
non-interacting part we choose to consider is as usual : 

^ m 

-:<™(^)^^V'L,™(a:):}. (26) 



Translation invariance along the y direction requires 
F{m' -\- s,m -\- s]n' -\' s,n + s) = F{m' , m; n', n). It is 
easy to write down RG equations for the F function to 
one-loop order, and we get : 



dtF{m' ,m;n' ,n) ^ > < F{rn\rn'';n'\n)F{rn'\rn;n\n'') ~ F{rn\rn'';n' ,n'')F{rn'' ,rn;n'' ,n)> . (27) 



rt\'\n" 



V 



A generalized Luttinger fixed point in real space is ob- 
tained for Fim! ,m]n' ,n) = Sm',m^n',nF(rn — n). This 
implies that the total particle number in each chain is 
conserved. In the vicinity of such a generalized Lut- 
tinger fixed point, we may write : F{m\m;n',n) = 
Sm' ,mSn' ,nF{m-n)+ f{m' , m; n', n). Here /(m', m; n', n) 
is assumed to be small in comparison to F{m — n). Lin- 
earizing the flow equations gives : 

dtf{m',m;n',n) = —B(m',m]n',n)f{m',m;n',n) , 

(28) 
with B(m' , m; n' , n) = F{m' — n') + F{m — n) — F{m' — 
n) — F(m — n'). The Luttinger fixed point characterized 
by the function F is locally stable if for any four-uple 
(to', to, n', n) we have 



B{m' , m\ n' , n) > . 



(29) 



Taking to,' = ?n + 1 and n' = n + 1 yields : 2F{m — 
n) - F{m - n + 1) - F{m - n - 1) > 0. But if 
m' = m -{- 1 and n' = n — 1, the same condition gives : 



2F(m - n + 1) - F{m - n + 2) - F{m - n) < 0. So it is 
clearly impossible to find a function F such that all the 
inequalities of Eq. (p9|) would be satisfied. In more physi- 
cal terms, real space Luttinger liquids are always destabi- 
lized at low energies by interaction processes which don't 
conserve the total particle number in each chain. At this 
point, we should emphasize that our analysis is strictly 
speaking valid only in the vicinity of the free electron 
fixed point, i.e. when F(m — n) is always small. Other- 
wise, the generalized Luttinger liquid would have to be 
treated non-perturbatively (using bosonization), which 
modifies Eq. (pS) for the anomalous dimension of per- 
turbing operators. Within this non-perturbative analy- 
sis, some recent workLj by Vishwanath and Carpentier 
suggests the possibility of stabilizing a real space Lut- 
tinger liquid. But this requires very special properties on 
the F function. The main difference between real space 
and fc-space Luttinger liquids comes from the number of 
free parameters which are required to describe the fixed 
point interaction. In the first case, the F function de- 
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pends on one variable, whereas in the second case, G* 
involves two variables, which leaves more possibilities to 
achieve local stability of the "tomographic" fixed points 
(at finite N). 

If we have N coupled chains with periodic boundary 
conditions in the transverse direction, we can go from the 
real space to the fc-space description by Fourier trans- 
form, according to : 



l^ik)^^=Y: / dxe<'^^+^'"^)^l^ix). (30) 



The real space Luttingcr fixed point with the F function 
yields the following interaction Hamiltonian, written in 
fc-space : 

^-t = 7^ E F{5)pn,i,5{<i)Pi^,J,-si-<i) ■ (31) 



LN 

I,J,S,q^O 

So, according to our general notation, we get Gs{I, J) ~ 

^{^) — X]m=o^-'^P(~*'^'''''^"^/-^)-^(''^)' '^hich depends 
only on the transferred momentum 5 along the flat Fermi 
surface. In Eq. (pi|), we drop the restriction that 1+6 and 
J — 6 should belong to the set {!,..., TV}, since we have 
the natural periodicity in fc-space : c]^ j(fc) = cj^ j i ^(fc)- 
To conclude both discussions, we would like to stress 
that each form of Luttinger liquid is very unlikely to ap- 
pear as the true low-energy flxed point of a real two- 
dimensional system. For an anisotropic metal with flat 
regions on the Fermi surface, the main difficulty arises 
in considering the limit where the number of patches be- 
comes infinite, in the case of non-singular microscopic 
interactions. For quasi one-dimensional conductors, we 
have shown that there are always some unstable direc- 
tions around the corresponding Luttinger liquid Hamil- 
tonian, even if we do not consider inter-chain hopping. 
In the remaining part of this article, we shall therefore 
consider only strong coupling regimes, but still within 
the perturbative RG approach. 



IV. STRONG COUPLING REGIME : THE RPA 

In this section we shall focus on the last two stages 
of a typical RG flow as shown on Figs. |l^ and |lj. The 
first task is to understand the plateau observed for in- 
termediate RG times. The comparison between relative 
contributions to the flow of Peierls and Cooper channels 
(see Fig. |5]) clearly shows that the former is dominant 
in this regime, with subleading corrections of order 1/A^ 
(see also Figs. |l^,|l^ and Table ||). 



A. RPA FIXED DIRECTIONS AND THEIR 
STABILITY 

It is first useful to present a complete description of 
the set of fixed directions within the RPA. In a sec- 



ond step, a local stability analysis around these direc- 
tions will be sketched. The RPA considered here consists 
simply in dropping the contributions from the Cooper 
channel in the RHS of the RG flow Eq. (^. As il- 
lustrated on Fig. pq, the notation is greatly simplified 
if we replace expressions like Gj-i{I,J + 9) by quan- 
tities of the form Fg{I,J), since 6 is the component of 
the transferred momentum along the flat Fermi surface, 
and is therefore conserved in the elementary particle-hole 
bubble. In the same spirit, we decompose the possi- 









K 



>Fs{I,J) 



J + t 



j + e 



FIG. 20: On the left, we show the interaction process with a 
notation adapted to the RPA. On the right, we give the sec- 
ond order contribution in the Peierls channel which involves 
Gk-i{I,K + 6)Gj-k{K,J -\- 6), and which is rewritten as 
FeiI,K)Fg{K,J) 

ble particle-hole states into singlet and triplet subspaces, 
which corresponding couplings are renormalized indepen- 
dently. More specifically, we set : 



rFrs(/,j)^G}_,(/,j 



3G}_,(/,J 



(32) 



1 F*"p(/, J) = G}_,(/, J + e)- G^j_j{i, J + e) 

With these definitions the RPA approximation reads : 
dtFo^{I,J)= J2 FS{I.K)FS{K,J). (33) 

1<K<N 

1 < K + e < N 

Here a denotes either the singlet or triplet channel. This 
has the generic form dtM = A/^, where M is a symmet- 
rical matrix of size Ng ~ N — 6, with real coefficients. 
It is straighforward to solve these simplified flow equa- 
tions for any initial condition M{t = 0) = Mg. There 
exists an orthogonal nratrix P and a diagonal matrix 
Ao such that Mq = PAqP^^. For any positive time, 
M{t) = PA{t)P-\ with [A{t)]u = [Ao],./(l - i[Ao]..)- 
In this situation, if some initial couplings [AqJm are posi- 
tive, the corresponding diagonal element diverges at time 
ti — [Aq]^^. Let us consider the set 1 of i's such 
that [Ao]m reaches its maximal positive value. For any 
ii,Z2 e I, [A{t)]i^i^ = [A{t)]i^i2. And it is clear that 
the ratio [A(t)]jj/[A(t)]ii vanishes in the limit where t 
approaches ti, provided j ^ T and i E T. This shows 
that the flow in M-space is appoaching a fixed direction 
along M* given by M* = PIjP^^. The matrix Ij is 
the diagonal matrix such that [Ii]m = 1 if i G X, and 
[Ix]m = if z ^ X. Let us denote by m the number of 
elements of X (we are assuming here that X is not empty, 
so m > 1 ; of course m < Ng). The fixed direction along 
M* can be interpreted as an orthogonal projector on a 
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m dimensional subspace of the Ng dimensional space on 
which M matrices operate. In this way, for a given m 
satisfying 1 < ?7i < Ng, we obtain a smooth manifold of 
fixed directions, which may be identified with the smooth 
manifold of m dimensional subspaces imbedded in a Ng 
dimensional vector space. This object is usually called 
a Grassmann manifold (G'm(K^'') being a standard no- 
tation) and its dimension is m{Ng — m). For a given 
fixed direction M* such that its largest positive eigen- 
value is m times degenerate, it is simple to characterize 
its attraction basin along the RPA coupling flow. Using 
a permutation of the columns of P, it is always possible 
to assume T — { 1 , . . . , m} . In the basis of M.^" defined 
by the columns of P, the set of matrices which flow at 
large enough times towards the direction dcflned by M* 
is given by : 



M = 



Im 

M 



(34) 



where Im is the m x m unit matrix, and M is any sym- 
metrical real matrix of size Ng — m whose eigenvalues 
are all strictly smaller than one. In this basis, M* cor- 
responds to Af = 0. So this attraction basin is an open 
set of matrices which is also a {Ng — m){Ng — m + l)/2 
dimensional manifold. The dimension of this attraction 
manifold is the largest when m = 1. This analysis shows 
we expect the Ng{Ng + l)/2 variables involved in the 
RPA flow equations in the vicinity of a ?7i-type flxed di- 
rection M* to fall in three classes (note that the notion 
of relevance or irrelevance refers to the behavior of the 
RPA flow in the space of directions and not to the global 
magnitude of the couplings) : 

- Irrelevant variables. We have {Ng — m){Ng — 171+ l)/2 
of them, and they correspond to changing the initial con- 
dition of the flow within the attraction basin of M* . 

- Marginal variables. They are associated to small 
changes of Mq within the Grassmann manifold Gm (M^" ) , 
therefore providing m{Ng — m) such variables. 

- Relevant variables. By a simple minded counting argu- 
ment, we expect their number to be m{ifn -t- 1)/2. 

This picture is conflrmed and refined by the following 
stability analysis. Let us consider a fixed direction M* = 
PljP^^ as above. For the sake of simplicity, we shall 
assume I — {1, . . . , m] (and thus Ij — I,„). This fixed 
direction is associated to the one parameter family of 
solutions M{t) = M*/{tc-t) of the equation dtM = M^. 
Now let us consider weakly perturbed solutions of the 
form M{t) = P{l,n/{tc -t) + N{t)}p-^. Linearizing the 
differential equation with respect to N{t), we obtain : 



dtN{t) 



IrnN{t) + N{t% 
tc-t 



(35) 



Solving this linearized equation gives : 

[N{t)],^ = [iV(0)],, (^)' if 1 < z,j < m, 

mt)h - [^(0)]., (ifri) if * < m and j>m + l 

or i > m + 1 and j < m, 



[ [N{t)h^ = [iV(0)]. 



if: 



l<hj <Ng. 



(36) 

The three cases correspond respectively to the relevant, 
marginal and irrelevant variables introduced before, jus- 
tifying the counting already given there. 

In the sequel, we shall often perform a local stabil- 
ity analysis around flxed directions of the form gi{t) — 
Ui/{tc — t) for flow equations of the form dtgi — AijkgjQk- 
This can be done by setting gi{t) = Ui/{tc ~ t) + Sgi{t), 
and by linearizing the evolution equation with respect to 
Soi. We obtain : 



dtSgi = f{t)B^kSgk, with 



f{t) = {U - t)-i 



(37) 

The linear stability is governed by the eigenvalues of the 
matrix B, which explicitly depends on the flxed direc- 
tion. Let us assume B can be put in diagonal form, with 
eigenvectors v^"\ such that Bv^°''^ — \aV^°'\ If Sg{t) is 
along v''°'\ that is 8g{t) — G{t)v^°'\ the linearized equa- 
tion implies that : g{t) = g{0)[tc/{tc ~ t)]^°. Setting 
e = 5(0), we have : 



9i{i) = f{t) { Ui + stc 



tc 



tr-t 



Aq-1 



,(a) 



+ 0{s') 



(38) 
The flxed direction along Ui will be stable if the term pro- 
portional to Ui grows faster than the one proportional to 

Vj , that is if all the eigenvalues Aq of B are smaller 
than 1. In the language already introduced, eigendirec- 
tions with Aq, > 1, resp. Aq = 1, resp. Aq < 1 will be 
referred to as relevant, resp. marginal, resp. irrelevant. 
Note that u is always an eigenvector of the B matrix 
with A = 2. This can be interpreted simply by a shift in 
the value of tr : 



Ui 



tc + Stc — t tc — t 



1- 



Str 



tc-t 



0{Stl 



(39) 



Applying this terminology, we find the A spectrum for 
these RPA fixed directions is very degenerate, with only 
three eigenvalues : A = 2 for relevant variables, A = 1 
for marginal ones, and A = for irrelevant directions (re- 
member Eq. (|3q)). Note that for a fixed direction with 
m = 1, we have only the trivial eigenvalue at A = 2, 
whereas the marginal A = 1 is Ng — 1 times degener- 
ate. We should mention that simiiap-stability analyses 
have been given by various authortSEj in the context of 
coupled chains of interacting fermions. 
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B. INFLUENCE OF THE COOPER TERMS ON 
THE RPA FLOW 

After this detailed description of the flow pattern in 
the RPA, let us describe the influence, on this flow, of 
a small perturbation in the Cooper channel. Intuitively, 
the Grassmann manifold Gm (K^" ) of fixed directions will 
evolve into an attracting manifold along which the flow in 
direction space is slow, by comparison to the variation of 
both irrelevant and relevant variables. In the remaining 
part of this section, we shall briefly discuss how we may 
construct this slow manifold using perturbation theory 
with respect to the Cooper channel. Then the dynamics 
along this manifold will be explicited. We shall consider 
the flow in direction space, of the form : 



dshi 



■^ijk f^j f^k 



7V2 



{Ajkihjhkhi)hi 



(40) 



For convenience we have set ||ft.|| = N instead of unity. 
This is because fixed directions in the RPA obtained from 
Ui = A^KJ^UjUk satisfy ||m|| = A^ as will be seen below. 

Now let us consider the flow with Aij]. = A]-], + eA]-],, 

where A'^^'^ refers to the Peierls channel and A^^^ to the 



Cooper channel. We are interested in the situation where 
e = 1, but it will be convenient to treat e as a free small 
parameter. Let us choose a fixed direction Ui — D{I, J) 
for the RPA flow, of the m = 1 type according to the 
above-given classification. The matrix D is then a pro- 
jector on a one-dimensional space, so there is a A^ di- 
mensional vector (j>{I) such that D{I, J) = (f){r)(f>{J) 

and J2i=i 't'il) — N- Here, for simplicity, we con- 
sider the case 9 = so Ng — N. We then have 
J2i jD{I,J) = N"^, so the norm of this vector D of 
couplings is equal to A^. 



To construct the point of the slow manifold which 
evolves smoothly from this fixed direction as s is 
switched on, we have to take a small perturbation svi — 
e{D{I,J),C{I, J)) which belongs to the irrelevant sub- 
space (i.e. the eigenspaceof i?(°^ with A = 0, according to 
the previous paragraph), incorporate Ui + evi in the RHS 
of Eq. (^) , and enforce that the projection of the result 
on the A = subspace vanishes, neglecting terms of order 
0(£^). We areJoUowing here the approach developped 
by N. FenicheEJ (see also the work by Van Kamper£3). 
We therefore obtain : 



J 



b{i, J) = ^HimJ) [^{if + HJf -wEk ^(k) 

, C(/,J) = -i0(/)'0(J)' 



(41) 



We may easily check that J2j D[I, J)4'{J) = 0, so £) be- A = 1) subspace. This procedure gives 
longs to the A = subspace. These expressions show 
that C(/, J) terms of relative strength 1/N are gener- 
ated, as well as a small difference between diagonal and 
off-diagonal terms of the D(/, J) matrix. These general 
trends will be confirmed by the explicit examples of fixed 
directions (including the Cooper terms to all orders in e) 
given in the following section. 



9,I?(/,J) = -— 0(/)0(J) 



^{ly + <t>{jy 



-j^Y.m' 



(42) 



It is easy to interpret this as a small evolution of 0(/), 
so that 



Now according to the intuitive picture already given, 
we expect that Cooper terms will induce a slow drift 
motion of the 4>{I) vector which doesn't depend on 
time if £ = 0. To compute this slow motion, we sim- 
ply have to consider the Cooper term for the unper- 
turbed fixed direction m^, namely to compute AiAujUk — 
■^{A-i!iUjUkUi)ui, and project this on the marginal (i.e. 



dsD{I, J) = dsCJ){I)cj){J) + Hl)ds^iJ), with (43) 



^^'^W - -^^-^W 



m'-^j:^iLf 



(44) 



Note that the norm of (^(/) is unchanged under this slow 
motion since we have J2i 4'{I)ds4>{I) = 0- 

Eq. ( |4^ ) has two interesting consequences. First, it is 
clear that fixed directions no longer form a continuous 
space, as was the case with flow equations involving only 
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a single channel. According to Eq. (p4[), fixed directions 
are obtained only if the p non-vanishing components of 
(j>{I) all have the same absolute value, equal to y^N/p. In 
the vicinity of such a fixed direction, it is easy to check 
that the tangent space to the slow manifold splits into 
p—1 eigenvalues A = 1 — Ae/p^ + O(e^) and N ~p eigen- 
values A = 1 -I- e/p^ + O(e^). The special case where 
p = N and 4>{I) — 1 independently of / is specially in- 
teresting, since it corresponds to an effective interaction 
with SU(A'^) symmetry. A detailed analysis will be given 
in the next section. We have represented what happens 
when iV = 2 on Fig. |^, in the (0(1), (/)(2)) plane. The 



order unity. With the normalization of this section, we 
have /ij(i) = iV5,(i)/||g(i)|| and dsjdt = \\g{t)\\lN. So 
the evolution of 119(^)11 is given by : 




FIG. 21: Flow pattern in the (0(1), 0(2)) space, in the Af = 2 
case. 

circle is of radius a/2, and is the continuous set of fixed 
point in the RPA. The small filled (resp. non-filled) cir- 
cles are the stable (resp. unstable) fixed points, when the 
Cooper channel is turned on. The arrows represent the 
direction of the flow. 

As a second important consequence of Eq. ([44[), it is 
now possible to estimate the duration of the intermedi- 
ate regime of the complete RG flow (see Figs. nG to 13). 



As already discussed earlier, the dominant couplings in 
this regime are of the -D(/, J) type. A very natural way 
to interpret the plateau regime seen on Figs. [lC| to p^ 
involves the slow motion just derived here. At the be- 
ginning of the plateau, -D(/, J) couplings are of the same 
order of magnitude, independently of J and J. But as the 
fictitious time s increases, there is a trend to eliminate 
couplings for which / or J is not one of the end points 
of the flat Fermi surface segments. In the language de- 
velopped in this section, the direction of the 0(/) vector 
evolves slowly in this regime under the residual influence 
of the Cooper channel. During this phase of the flow, 
the magnitude of the coupling vector increases tremen- 
dously. Let us estimate this semi-quantitatively. From 
Eq. (Q), we see that we have to wait for a s interval of 
order A^^ to get a change in a given component 0(/) of 



3sln||g|| == -TT^Aijkhihjhk ■ 



(45) 



To get an estimate, we approximate the RHS by its value 
obtained if e = 0, namely keeping only the Peierls chan- 
nel, and assuming D{I, J) = (j){I)cj){J), C{I, J) — and 
J2 4>{I) = N- With these approximations, ds In ||g|| ~ 1, 
or |jg(s)|| = ||gi(0)|| exp(s). Since s increases by an 
amount of order A''^, we see that the magnitude of the 
coupling vector is multiplied by a factor of the form 
exp(7A'^^), where 7 is a positive number which is expected 
not to depend on N. Unfortunately, this restricts tremen- 
dously the possibility to observe the flnal stage of the flow 
in physical systems. Indeed, our whole approach relies on 
a one-loop approximation, and therefore requires ||g||/A^ 
to remain small. This criterion breaks down as s goes to 
inflnity. But even to observe the complete plateau regime 
in a physical system would require extremely weak ini- 
tial couphngs, of the order of exp(— 7^^), which doesn't 
make much physical sense, since N is in principle infinite. 
This shows that for all practical purposes, the system is 
in many respects "frozen" at the entrance of the plateau 
regime, therefore justifying the description already given 
by Zheleznyak et al. Note however that we cannot rule 
out at this stage the possibility of a cross-over from a 
spin density wave state to a d-wave state superconductor 
at lower energies which would be induced by higher order 
terms in the RG fiow, beyond one-loop. But we don't see 
any strong argument to support this scenario. Any at- 
tempt to investigate higher order corrections is severely 
limited by the divergency of the couplings in the lowest 
order approximation. 

In the following sections, we shall study in great detail 
some fixed directions where (j){I) is non-vanishing for p 
values of /, focusing on extreme cases p = N and p = 2. 
p = 2 corresponds to the final regime of the one-loop 
fiow, whereas p = A^ is a particularly interesting fixed 
direction, since the flow passes in a rather close neighbor- 
hood of this point at the beginning of the intermediate 
regime. Its investigation is greatly simplified by the fact 
it exhibits a global SU(A^) symmetry for the effective 
couplings. 



V. STUDY OF FIXED DIRECTIONS THAT 
ARE HIGHLY SYMMETRICAL 

A. DETERMINATION OF THE SYMMETRIC 
FIXED DIRECTIONS 

We now aim at studying in greater detail the interme- 
diate regime, for which the couplings show an approx- 
imate SU(A'^) symmetry. As usual, the use of a sym- 
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metry greatly simplifies the calculation, so that we will 
restrict ourselves to the case where the couplings are ex- 
actly SU(A^) symmetric. 

Fixed directions are found when solving the equation 
Ui — AijkUjUk, where we have used the general notation. 

I 



In the particular case we are interested in, we simply have 
to consider Eq. (16), and replace dtG in the LHS with G. 
Considering only the princip al co uplings of the type A, 
C and D, as discussed in Sec. [IC, and setting all others 
to zero, we get the following system : 



NA"" = AA""^ + 2{N - 1)(D"^ + D^Z?") 

ND"" = {N- 2){D''^ + 31?"^) + 2{A' - C')D'' + 6(A" - C'')D'' 
^ ND"" = 2{N - 2){D''^ + D^'D'') + 2{A^ ~ C'')D'' + 2(A'= - C'')D'' + 4(74" + C")D'' 



(46) 



r 



The SU(iV) symmetry is imposed by the relations : in Eq. (ES), in the charge and spin sector, is then of the 
^c,s _ ^c,s j_ 2^c,s jQ(jgg(;i^ tj^g interaction Hamiltonian following schematic form : 

I 



q^O I,k 



.J,k' 

-^ E [E4„/(fc+9K,/(fc')][E</(^'-9)^R^Kfc)' 



q^O,k,k' J 



(47) 



TABLE II: Values of the couplings for the SU(iV) fixed direc- 
tions. 





Spinless 


(+,+) 


{+,-) 


A" 


1-1/iV 


(1 - l/iV)/4 


3/4[l - 1/N - 2/N'' 

+0{l/N')] 


c 


-1/iV 


-l/(4iV) 


-3/{4N)+0{l/N-') 


D" 


1 


1/4 


3/4[l - 2/N'^ + Oil/N")] 


A' 





1/4 


-1/4[1 -2/N + 2/N^ 
+0{1/N')] 


C= 








1/{2N)+0{1/N-') 


D= 





1/4 


-l/Al + 2/N'^ +0(1/N'') 



TABLE III: Values of the couplings for the SU(2) fixed direc- 
tions ("hot spots" case) . 





Spinless 


(+,+) 


{+,-) 


A^/N 


1/4 


1/16 


1/16 


cyN 


-1/4 


-1/16 


-1/16 


oyN 


1/2 


1/8 


1/8 


A^'/N 





1/8 





C/N 








1/8 


oyN 





1/8 


-1/8 



The spinless case is obtained when setting all spin cou- 
plings to zero. A little bit of algebra leads to the results 
listed in Table O In the case of spin 1/2 electrons, we 
find two fixed directions, named (-I-, -I-) and (-I-, — ), where 
the signs are those of {D'^, D^). In the spinless and (-I-, -I-) 
cases, the results are exact, whereas in the (-1-,—) case, 
we gave a 1/N expansion of the results valid up to order 
3. In fact the ZJ's are solution of a polynomial equation 
of order 4 and thus have complicated exact expressions. 
The (-I-, -I-) fixed direction corresponds to a charge den- 
sity wave, whereas the (+, — ) fixed direction corresponds 
to a spin density wave. This will be shown below, when 
we compute the response function in these cases. 



Another simple and important case is the fixed direc- 
tion that involves only the 2 patches at the border of 
the Ferm i surface, that we have called "hot spots" in 
Sec. [IC. We have indeed seen that this fixed direction 



was reached in the final part of the flow, after the flow has 
spent quite a long time near the (+,— ) fixed direction. 
We have looked for such "hot spots" fixed directions, and 
we have found three of them, respecting a SU(2) symme- 
try. The resulting couplings (divided by N) are shown in 
Table [IF Notice that most of these results can simply 
be obtained by setting A^ = 2 in Table p. The only ex- 
ception is in the (-I-, — ) case, where the whole 1/A^ series 
would have to be summed. 
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B. RESPONSE FUNCTIONS FOR THESE 
FIXED DIRECTIONS 

1. Notations and general setting 

The formalism of the renorniahzation of the response 
fuBctions is well known (see for example Zheleznyak et 
alll3), but let us briefly summarize it and fix our nota- 
tions. We shall here use the RG to follow the evolution 
of the response functions with the decreasing of a typical 
energy scale i'. The high energy cut-off Aq is fixed, and 
here we set t — ln(Ao/z^), so that t increases when i/ is 
lowered. The RG equations for the couplings are just the 
same as before (see Eq. (E&)). In this case, t/ is the energy 
scale at which the interaction process takes place. We 
now consider the response functions. We introduce the 



triangular vertices and the susceptibilities (sec Fig. 22). 
From the point of view of a renormalizable field theory, 
these are the only composite ooerators that need being 
renormalized (see for instanceE3). More physically, the 
triangular vertices describe the interaction between the 
fermions and an external field, whereas the susceptibili- 
ties give the linear response of the system when a small 
external field is applied. When a susceptibility diverges, 
the system gets spontaneously ordered, even in a vanish- 
ing external field. The common wisdom, when several 
susceptibilities diverge, is to assume that the most di- 
vergent susceptibility determines what phase transition 
is observed. 



J + 5 

--< — 
J 

J 





'\jyrvf\JD. 



n 



■rrii^J) 



4 




Kf(/,J) 



FIG. 22: Representation of the triangular vertices and sus- 
ceptibilities in the particle-hole and particle-particle channels. 
We have used the notation /^ = A^^-l-l — 7-|-5 and similarly for 
J'g. I's^o would be equal to —7 if the patch numbers ran from 
— 7V/2 to A''/2 instead of 1 to A^. For the sake of simplicity, 
we have not mentioned whether the singlet or triplet channel 
is involved, and we have only indicated the momenta parallel 
to the Fermi surface (i.e. the patch numbers), but not the 
perpendicular ones, nor the energies. 



As is usual, we introduce different response functions, 
corresponding to different physical situations. These 
can be either of the particle-particle (Cooper) or of the 
particle-hole (Peierls) type, i.e. they concern either the 
superconducting instability or the density wave instabil- 
ity. In our notation, we will associate pp and p*^ super- 
scripts. This is shown on Fig. |23. Notice on this figure 



TABLE IV: Dictionary between the charge-spin and 
singlet-triplet ways of separating the couplings. 



the 





particle-hole 


particle-particle 


singlet 


G" + 3G= 


G" - 3G= 


triplet 


G'-G' 


G' + G'' 



that we have only indicated the momenta parallel to the 
Fermi surface, that is the patch numbers. None of these 
momenta is summed over, contrarily to the momenta per- 
pendicular to the Fermi surface, and the energies that 
appear in the loops. For the external legs, we choose the 
configuration that gives the biggest divergence. We thus 
choose all external fermion momenta on the Fermi sur- 
face, and the momentum entering the vertices are chosen 
equal to 2kp (resp. 0) in the particle-hole (resp. particle- 
particle) channel. The energy entering the vertices is 
chosen to be v as discussed above. 

We insist it is not sufficient to tell what comes in and 
out in the triangular vertices (the patches involving J, 
J +S or J'g). It is necessary to keep track of what happens 
at the vertex pictured by a circle, i.e. of the non-diagonal 
character of the triangular vertices. If we consider the 
patch number as a tag for different fermion species, this 
phenomenon isJcnown as operator mixing in the field the- 
ory literature.E3 Notice also that since our patch numbers 
run from 1 to A^, it is /^ = A^ + 1 — / + (5 that appears, 
instead of — / -|- (5 as would be the case if the patch were 
numbered from —N/2 to N/2 (in the even A^ case). 



We furthermore distinguish between the singlet and 
the triplet channels, respectively giving rise to singlet and 
triplet superconductivity in the particle-particle channel, 
and to charge and spin density wave in the particle-hole 
channel. We will add a s or t subscript in our notation 
for the response functions to indicate what channel is 
under consideration. For example, the spin density wave 
susceptibility is denoted by TZ^ ^{I, J). From the spin ro- 
tation invariance of the interaction, it is clear that singlet 
and triplet channels do not mix. It is thus interesting to 
separate the interaction process in singlet and in triplet 
channels, instead of the charge and spin couplings. This 
is done differently in the particle-hole and the particle- 
particle cases, because the rotation group doesn't act in 
the same way in both cases (the representations are differ- 
ent). The resulting couplings in the four distinct channels 
are given in Table f^. 

We now come to the RG equations for the response 
functions. The different Feynman diagrams that con- 
tribute to the renormalization of the different response 
functions are shown on Fig. E3. It can be shown that the 
most divergent responses in each channel always appear 
ior 5 — so that we will assume (5 = in all the follow- 
ing, and we will not write the 6 subscript anymore. The 
RG equations in the singlet and triplet channels are the 
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K + S J + S 




-X/^iX/XjO' 



^\jyf\j^^. 




;a/\Vrv/^ 



;a/\S^v/^ 



FIG. 23: Feynman diagrams that contribute to the renormal- 
ization of the triangular vertices (left) and the susceptibilities 
(right) in the particle-hole channel (top) and particle-particle 
channel (bottom). The K variable is summed over. 



TABLE V: Biggest exponents of the susceptibilities, for the 
four fixed directions involving spinning electrons. 





SU(iV) 


SU(iV) 
(+r) 


SU(2) 


SU(2) 


ph, s 


f-l/(2iV^) 




7/8 




ph, T 




l-9/{2N^) + 0{l/N'') 






pp, S 








7/8 


PP, T 











can easily be cast in a matrix form : 
dtTPP = TPPTWPP 

9J7^ph = r^^ trp^ 
a7^pp = Tpp 'Tpp 



(49) 



same and read (we drop the s and t subscripts here) : 
f dtTP^il, J) = ^ Ea' ^p'^U, K)Gj-k{K, J) 

dtW'-il, J) = Y.K '^'"'{1, K)TP''{J, K) 

y dtnppii,j)^j:j^TPPii,K)TPP{j,K) 

(48) 
where K' — N +1 — K. It should be noticed that in fact, 
there should be two different triangular vertices in each 
channel. Let us consider the particle-particle (and for ex- 
ample singlet) channel. The two external fermions can go 
towards the vertex or go away from the vertex. It is how- 
ever quite easy to see that those two vertices are related 
through complex conjugation. But here our vertices are 
all real. Indeed they are real at i = (i.e. v = Aq), 
their value being their non-interacting one, that is 1 for 
the diagonal elements and for the others. Furthermore, 
from the renormalization equations that do not involve 
any complex numbers, they have to remain real along the 
flow. Let us also remark that the initial values of the TZ^s 
are all zero, which is their non-interacting value, when 
evaluated at i' = Aq. In Eq. ([48|), the G's are of course 
channel dependent, and are found thanks to Table IV. 



2. Analytical results 

The previous equations are directly usable in a numer- 
ical simulation, once the t variable is changed for the s 
variable, which is simply done by dividing these equations 
by the norm A*" of the coupling vector. Before showing 
numerical results, it is interesting to give analytical re- 
sults that are available for the fixed directions we found in 
Sec. V A. Let us briefly explain how this works. Eq. (E3) 



where *T is the transpose of T. The definitions of the 
Al's are clear, and they are symmetric thanks to the 
particle- hole symmetry of the couplings (coming from the 
hermiticity of the Hamiltonian) : Gs{I,J) ~ G-s{J,I)- 
For a fixed direction, we know that Ai — U/{tc — t), 
where lA is a t independent matrix that defines the fixed 
direction (it is formed with the Mi's). We diagonalize this 
matrix, setting U = P'DV~^, where V is diagonal, and 
V is an orthogonal matrix. Denoting T = V^^TV, and 
similarly for the TZ's, we get : 



StTP^ ^ -i-T Phpph 



dtTPP = 
9t7^P'^ = 
dtTZPP = 



_!_7'PP2?PP 

'T'ph t'T'ph 
-j-pp t-j-pp 



(50) 



From these equations and the initial conditions T{t = 
0) = I and TZ{t = 0) = 0, it is evident that the T's and the 
7?.'s are diagonal at any time. One of the consequences 
is that T = *T. Let us suppose that [Dju — 6i ^ 1/2. 
It is easy to check that one gets : 






tr 



tc-t 

tc 

25, - 1 



, and 



tr -t 



2(5; -1 



(51) 
(52) 



If 5i = 1/2, we get a logarithm instead of power functions 
for the susceptibilities (which is the case in the purely ID 
case). In order to get the critical exponents 2^^ — 1, we 
thus only have to diagonalize a matrix. We have com- 
puted these exponents, and in Table M we give the biggest 
exponent of the four channels, for the four fixed directions 
we have previously found. All the other exponents can 
be shown to be equal to -1, up to some 1/A^ corrections. 



The results in the spinless case are given in Table VI 



We also find the other exponents are nearly equal to -1 
in the SU(iV) spinless case. We thus see that in the 
SU(A^) case, the (+,+) fixed direction is of the charge 
density wave type, whereas the (+,-) fixed direction is of 
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TABLE VI: Biggest exponents of the susceptibilities, for the 
two fixed directions involving spinless electrons. 





SU(iV) 


SU(2) 


ph 


I - 2/iV^ 


1/2 


pp 




1/2 



the spin density wave type. No superconducting insta- 
bility arises in this case. In the SU(2) case, i.e. the "hot 
spots" regime, the (+,+) fixed direction is still of the 
charge density wave type, but the (+,-) fixed direction is 
of the singlet superconductivity type. Let us emphasize 
that even if the exponent 7/8 in the SU(2), particle-hole, 
singlet, (-l-,-l-) case, can be found in taking iV = 2 in 
1 - 1/{2N'^), the iV = 2 case is different from the iV 7^ 2 
case, as can be seen from the (-1-,-) fixed direction where 
the dominant instability is not of the same type. In the 
spinless case, the SU(iV) fixed direction gives rise to a 
charge density wave, and the SU(2) fixed direction ex- 
hibits a perfect competition between charge density wave 
and superconductivity. 

In all three SU(A^) fixed directions, the most diver- 
gent susceptibility is obtained by summing all the sus- 
ceptibilities : TZ* = J2i j'^i^' ■-^)- This corresponds to a 
uniform sign of the particle- hole or particle-particle wave 
function. This is also true in the SU(2) case, except for 
the (-1-,-) fixed direction, where the most divergent sus- 
ceptibility is found to be : 7^* = 7^(l, 1) - 7^(l,iV) - 
TZ{N, 1) + Tl{N, N). This can be checked to correspond 
to a "d-wave" type particle-particle wave function. 

The d-wave superconducting correlations found for the 
SU(2) (-f-,-) fixed direction looks similar to the CISO 
phase described by various authors for doped 2-leg Hub- 
bard laddfl rs . Sin ce ladder systems have been extensively 
studiedJ^'r^lrHr'l it seems useful to explain the connec- 
tion between our systems and 2-leg ladders. In a 2- 
leg ladder, the single-particle spectrum is composed of 
two energy bands with the following dispersion relation : 
£a{k) — —2tcos{k) — ij^cos(a). Here k is the momen- 
tum along the ladder, and a £ {0, tt} is the transverse 
momentum. If t± is not too large and the filling factor 
not too large nor too small, both bands cross the Fermi 
energy, at wave- vectors fcp.a- If t±_ is positive, we have 
kF,-K < kpfi- From the fact that /cp.ir 7^ ^f,Oj the D 
processes become irrelevant at low energy, as illustrated 
on Fig. ^(a). On the other hand, this system allows 
for transverse Umklapp processes, which are depicted on 
Fig. |2^(6), and these will be denoted here by F'^'^. 



kr.Ti 



kp,; 



D' 



kp. 



r '°'w- 



fC,S 



k 



F,0 



FIG. 24: (a) The D'^'^ processes are no longer available at low 
energy on a 2-leg ladder since fcp.o 7^ fep^jr. (b) But transverse 
umklapps are now possible ; they are denoted by F'^'". 



Of course, ^'=•''(0), A^'^tt), C"='"(0,7r) = C"='"(7r,0) = 
(^c.s pj-Qcesses remain available. To simplify the com- 
parison with our model, we shall assume the 2-leg lad- 
der is not half-filled so that longitudinal Umklapp pro- 
cesses are irrelevant. We shall also neglect the differ- 
ence between the two Fermi velocities, and assume that 
^'"■"(0) = ^'=■''(71) = A'''''. Clearly, our model of cou- 
pled chains in fc-space in the hot spot regime, where only 
chains labelled by 1 and N are involved in the dominant 
scattering processes is not directly equivalent to an effec- 
tive model for a doped 2-leg ladder. The D'^'^ couplings 
in the former case are replaced by the F"^'^ ones in the 
latter. However, it is possible to connect both situations 
by a unitary transformation U. We shall define U as 
follows : 

UcnjAkp-' = CR,/,^(fc), (/ G {1, A^}) (53) 
Ucui,^^i^k)U-' - ficlijk - 2fcF). (54) 

Here, we have assumed that the spin index fi takes two 
possible values ±1. Using the definitions for the charge 
and spin currents (Eqs. |lj and |l^), it is clear they are 
unchanged by U on the right-moving side, but they are 
transformed on the left-moving one according to : 



usijAq)u- 



-Phj+s-s{q), 
Sh.i+s.-sil)- 



(55) 
(56) 



Furthermore, because of particle-hole symmetry, U com- 
mutes with the kinetic energy. Starting from a model 
with couplings A'^''^, C^'^, and D'^'^, and after applying 
the unitary transformation U, we get a model with cou- 
plings A'^'^, C"^''', and F"^'^, given by : 

A'' = -A''; & = -C; F" = -D" and 

i" = A"; C" = C"; #" = L>" (57) 

Our (+,-) SU(2) fixed direction can be reached directly 
from the particular choices of initial couplings such that 
A'^,^{1) = A'^'^'iN) = G"^■^ C"='''(1,A^) = C"^'"(iV,l) = 
C"'^, and i:>'='"(l,iV) = D'='''{N,1) = G"=^% aU other bare 
couplings being set to zero. After the U transformation, 
we get a model equivalent to a doped 2-leg ladder. In the 
g-ology notation used by Schulz,t3, the bare couplings are 
.91 = -2G" and 32 = -G= - G^ The {+,-) SU(2) fixed 
direction is the stable low-energy attractor when C^ > 
and G^ < 0, which translates into gi > max(0, 2na) in 
the 2-leg ladder terminology. According to Schulz,E3 the 
2-leg ladder is in the staggered flux phase regime (see 
Fig. 1 in Schulz's paperc3). This is in perfect agreement 
with our d-wave superconductivity dominant correlation 
for the (-1-,-) SU(2) fixed direction. Indeed, the d-wave 
order parameter reads : 

k 

+ 4l^N,li>^F + k)cl,^^{-kp-k). (59) 
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The unitary transformation U turns this into 

y^ c|^,i,o-(fcF + fc)cL,A',<T(-fcF + k) 



k,cr 



-R,N,iT 



(fcp + fc)cL,l,(T(-^F + k). 



(60) 



To translate this in the language of a 2-leg ladder we 
assume that chain 1 corresponds to a = and chain 
N to a = TT. Shifting a little bit the fcp values, which is 
harmless since only the marginal couplings at zero energy 
have been considered, we get for the order parameter : 



y^ CRfi.aikF,0 + fc)cL,7r,o-(-A:F,Tr + k) 
fc,cr 

~'^R,7r,o-(^F,7r + fc)cL,0,o- (-A:f,0 + k) . 



(61) 



In real space, this reads 



E' 



,i(feF,o+feF.,r)a;r /t 



[4{x)M^) - 4ix)M^)l (62) 



where ipj ^{x) creates an electron with spin a along the 
rung X of the ladder, and on chain number i. Clearly, 
this corresponds to a staggered current pattern. 

This phase has many interesting properties. It has 
only one gapless mode, which is the even charge riHode. 
All other modes are gapped. As shown by Schulzjld this 
system exhibits a large SO (6) symmetry, which contains 
our SU(2)xSU(2). But a detailed analysis of the energy 
spectrum of this model is beyond the scope of this article. 



3. Numerical results 

In this section we will show some numerical flows of 
the various response functions TZ, for the spinful case, 
and for initial conditions such that the couplings get close 
to the SU(7V) (+,-) fixed direction, in the intermediate 
regime. We consider a number of patches N — 8 and 
not A'' = 16 as in Figs. 10 to O. The reason is only to 
shorten the computation time. Indeed, the final growth 
of the "d-wave" superconducting response function does 
not occur immediately after the transition from SU(A) 
to SU(2), but only after a delay time that is growing 
with A^. The flows are shown on Fig. ^. In Fig. |25|(a) 
we show the evolution of the D'^ couplings, exactly as 
in Fig. nq, so that one can read the times at which the 
intermediate approximate SU(A'^) regime is reached, and 
at which the transition to the SU(2) fixed direction takes 
place. The flows of the biggest response function in each 
of the four channels (particle-particle/particle-hole and 
singlet/triplet) are represented in Fig. ^(t*). We clearly 
see that during the intermediate regime, the biggest re- 
sponse fuction is in the particle-hole, triplet channel, in 
other word a spin density wave. We have checked that it 
corresponds to an approximately uniform "s-wave" type 
response function. It would be exactly uniform only if the 




FIG. 25: Flows of response functions for N — 8 and initial 
condition C^ = —G^ = 0.1. (a) : flow of couplings D'^ as in 
Fig. hol (fe) : flow of the biggest response function in each 
channel, (c) : zoom of (6) for small time s. (d) : competition 
at small time s, in the particle-particle, singlet channel, of the 
"d-wave" and "s-wave" types of response functions. 



SU(A^) symmetry was exactly realized. After the transi- 
tion to SU(2), this response function saturates while the 
particle-particle, singlet response function grows and be- 
comes bigger than the former. We have also checked that 
it corresponds to a "d-wave" type response function. We 
have shown in Fig. |2^(c) a zoom of Fig. [25|(5) so that one 
can see the initial competition between the four chan- 
nels. Finally, Fig. E3((i) shows the competition between 
two superconducting response functions. The first one 
(represented with circles) is of the "s-wave" type and is 
the biggest in the beginning of the flow, while the sec- 
ond one (represented with squares) is the "d-wave" type 
response function we have previously already mentioned, 
that becomes the biggest after a while. Even if those two 
response functions are much smaller than the spin density 
wave one, it is interesting to see how in their competi- 
tion, the "hot spots" already manifest themselves in the 
intermediate regime. 



C. PHYSICAL PICTURE FOR THE FIXED 
DIRECTIONS WITH SU(7V) SYMMETRY 



In this section, we shall attempt to discuss some 
properties of these fixed directions which are not easily 
grasped within the perturbative RG approach followed 
in most of this work. For instance, it has been known for 
several decades that many electronic systems exhibit spin 
and charge decoupling in one dimension. In the present 
case, the transverse momentum is an additional degree 
of freedom which we shall denote by "color" in this dis- 
cussion. A natural question is therefore : do we have 
decoupled charge, spin and color excitations in these sys- 
tems ? It is useful in this respect to introduce charge. 
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spin and color currents, which may be defined as 



PR (9) 



E 



kJ,X ^R,I,X 
k,I,\,p,^RJ,\\ 



(fc + q)cYi.i.\{k) for charge 
(fc + q)al CYij^f_,{k) for spin 



A ^R,/,A 



(fc + (7)cR.j_A(fc) for color 



(63) 



As it is well known, it is possible to express the kinetic 
energy part Hq of the system Hamiltoniaji as a sum of 
three quadratic terms in these currents C^ Hq — ifo.R + 
i/o.L, with : 



Ho.] 



nv-p 
MN 



XI • Pr(9)Pr(-9) 



■KV-p 

' M + N 



E 

qJJ 



TTVp 



2{M + N) 



E ■■ SRiQ)Sni-q) 



1 



Ji-'{q)Ji^'i-q):-j^:Ji'iq)Ji-'i~q) 



(64) 



r 



For the right moving branch, normal-ordered products 
are defined according to : A{q)B{q') := A{q)B{q') if 
q > q' and B{q')A{q) if g < q' . Note that similar ex- 
pressions may be written for the left moving branch. We 
have introduced M = 2 for SU(2) spins, to illustrate 
what would happen for SU(A/) spins. It is important to 



note that each term in i?o.R commutes with the other two 
terms. This property is often stated as charge-spin-color 
decoupling. Now, the most general interaction which 
preserves SU(2)xSU(A^) symmetry and this decoupling, 
may be written as : 



Hint — ffch. 



argo / , 



pi.{~q)pK{q) 



t/spm / ^ 
q,a 



s't{-q)sUq) -\- 9coio. 



q,I,J 



^l\-<DJi\<D-j^JU{~<i)J^\<D 



(65) 



r 



These interaction terms may be cast into the 
parametrization used in this section, and the precise 
mapping reads : 



^ ^ — .^charge ^ ^color 

D'^ — D^ = gcolor 

C^ — n ■ 
^^ — yspin 



Sspin 



(66) 



So the fixed directions with SU(2) x SU(iV) symmetry dis- 
cussed here cannot be written in this charge-spin-color 
decoupled form, since D'^ 7^ for both (+,+) and (+,-) 
fixed directions, as seen on Table ||. We do have a weaker 



form of decoupling, where the charge bosons prl {q) de- 
couple and remain in a gapless Luttinger liquid state. 
But spin and color sectors are strongly mixed. Actually 
the (+,+) direction exhibits a larger SU(27V) symmetry 
than the expected SU(2)xSU(A^). To show this, let us 
define generalized SU(2A^) currents by : 

Elt'''^{q) = Y. clj,^{k + q)cnjAk) , (67) 



and let us consider the U(l)xSU(2A^) invariant interac- 
tion : 



iJint = gcy..r,oJ2pL{-q)pRiq)+g ^ \ eI''^'-^ i-q)E^^''' '\q) :^E['''''\-q)Ei^^''^^iq)\ 



(68) 



r 



It is easy to show that such an interaction Hamiltonian couplings : 
corresponds to the following values of the charge and spin 



< 



^ ^ — .ychargc 27V ' 

D" - Z?" = 



(69) 



26 



From Table |l|, it is clear that the (+,+) direction can 
be written this way, with ijchargo — and g = 1/2. 
This model is simply the SU(27V) Chiral-invariant Gross- 
Neveu modeJ, which was diagonalized by Andrei and 
Lowensteinca using the Bethe Ansatz for many-body 
wavefunctions. From their results, a rather simple pic- 
ture emerges, involving a massive SU(2iV) "spin" sec- 
tor, completely decoupled from a single massless bosonic 
mode (the charge mode already discussed). In the TV -^ 
oo limit, this physical picture evolves smoothly towards a 
mean-field description, characterized by the appearance 
of a non-vanishing expectation value for the charge den- 
sity wave operator ^^ c^^ j -^{k + A:F)cL,/,A(fc — ^f)- This 
mean- field shows itself in the behavior of the particle-hole 
response function in the singlet sector, which diverges 
near the instability with an exponent equal to 1— 1/(27V)^ 
(see Table 0). As iV goes to infinity, the mean-field (or 
RPA) value 1 is indeed recovered. For large but finite N, 
long wavelength fiuctuations of the order parameter pre- 
vent a true spontaneous symmetry breaking. However, 
as discussed long ago by several authors ,e3lj some qual- 
itative conclusions of the mean-field analysis are correct, 
as for example the prediction of a non-vanishing energy 
gap in the single-electron spectral function. Note that 
the original electrons are no longer the fundamental ex- 
citations, and they split into a massless charge part and 
a massive neutral spinful part, so their spectral function 
does not exhibit a sharp quasi-particle pole. 

Let us now focus on the (-1-,-) fixed directions. As 
hinted by Eq. |47|, the corresponding interaction Hamilto- 
nian is of the form : 

Hint = g dx : Vr / A i^)<^\x' V'L,/. A' (a;) 
Jo ' ' 

^i^ljJ^Kf.'^RJ.^^'i^) ■■ ■ (70) 

In a mean-field picture, the corresponding order parame- 
ter is £_°-{x) =< '0R„7,A(^)^AA'V'L,/,A'(a;) >. ^{x) is a vec- 
tor with three complex coordinates. It transforms under 
the SO (3) rotation group corresponding to global rota- 
tions in spin space. Furthermore, under a chiral trans- 
formation which changes i/'r / a i^^*^ exp(i0/2)'0j^ j ^ 

and i/'l / A ^^^'^ 6'^p(~*^/2)V'l / A' ^(^) ^^ replaced by 
exp(z6')^(a;). Note that a global chiral transformation 
corresponds roughly to a uniform translation of the spin 
density wave condensate. So the mean-field approxima- 
tion breaks two independent continuous symmetries of 
the underlying Hamiltonian. It is therefore natural to 
derive an effective action for the long wavelength fluc- 
tuations of the order parameter ^{x). This problem has 
been addressed by several authors in the context of dy- 
namical properties of spin density waves. For a recent 
source, with several references to earlier wotk, see for in- 
stance the article by Sengupta and Dupuis.E^ In order to 
keep the presentation simple, we shall simply quote the 
results obtained in the long wavelength approximation, 
keeping only massless modes. It turns out that these 
modes correspond to order parameter configurations of 



the form ^{x) = pexp{i9(x))'n{x), where p is a positive 
number found by solving mean-field equations, 9{x) is 
the "phason" field, and n{x) is a real unit vector, giving 
the direction of the staggered magnetization. Integrating 
out fermions provides the following effective action : 

S^ff{0, n)^^ I AxAt {d^.ed^'e + d^.nd^'n) (71) 

-I- less relevant terms. 

As expected, the phason field remains gapless, and its 
fluctuations preserve the chiral symmetry of the system 
at low energy. Now, the magnetic local order parameter 
will also remain in a symmetric phase with respect to 
spin rotations. However, by contrast to the phason, it 
will develop a spectral gap M ~ ?Tiexp(— A^), where wis. 
the fermion mass gap (see for instance Polyakov's book,c£l 
chap. 2). M is therefore a spin gap. Its appearance can- 
not be detected in the simple perturbative approach in 
terms of the original fermi ons. Indeed, the correlation 
functions computed in Sec. VB above were diverging at 
the typical energy scale m at which particle-hole bound 
states occur. Clearly, our one-loop fermionic RG misses 
long wavelength fluctuations of these composite bosons 
completely. Coming to real systems, we see that this 
spin gap is expected to vanish quickly in the N —f oo 
limit, which again has to be taken in order to describe 
two-dimensional systems with a continuous Fermi sur- 
face (this exponential variation of the spin gap with N 
is very similar to the behavior found for an even num- 
ber of weakly coupled antiferromagnetic,spin chains; see 
for instance the work by ChakravartycZl) . To conclude 
this analysis, we emphasize that most physical proper- 
ties of the (-1-,-) fixed directions are well captured, as N 
becomes infinite, by a standard mean-field treatment of 
the Peierls instability. 



VI. STABILITY ANALYSIS OF THE 
SYMMETRICAL FIXED DIRECTIONS 

A. LINEAR STABILITY ANALYSIS 

In this section, we shall study the stability of the 
SU(7V) and SU(2) fixed directions studied in the previ- 
ous section. We have already studied stability questions 

This 



within the general setting of the RPA in Sec. IV 



stability analysis was performed around the fixed direc- 
tions obtained with a small contribution of the Cooper 
channel (remember we introduced a small e as a way to 
tune the influence of the Cooper channel with respect to 
the Peierls one). Here we directly study the stability of 
the SU(A^) and SU(2) fixed directions, for which no as- 
sumption about the smallness of the Cooper channel is 
made. The very possibility to compute the eigenvalues in 
a linear stability analysis in these cases comes from the 
gigantic simplifications gained from the symmetries. We 
recall that with the general notations (see Sec. IV), the 
eigenvalues are those of the matrix B (Eq. p7|). 
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We shall briefly explain how to compute the eigenval- 
ues and their degeneracies, in the SU(A^) spinless case 
only. The other cases are either much simpler (SU(2)), 
or the algebra is much heavier but the ideas just as sim- 
ple (SU(A^) with spin). The fixed direction only involves 
the principal couplings A, C and D. It can be shown 
that a perturbation in the non-principal couplings is ir- 
relevant, so that we will restrict ourselves to the principal 
couplings. We set D{I, J) = D/{tc-t) + SD{I, J), where 
D = 1 is the value of the D-couplings on the fixed direc- 
tion (see Table H). Doing the same with the A's and C's, 
linearizing Eq. ( |l7| ) (it would be Eq. ( [l^ ) in the case with 
spin) , and writing down the eigenvalue problem yields : 

AM(/) = ^Sil) 

X5C{I,J) = -^SD{I,J) 

X6D{I,J) = ^[S{I) + S{J) 

+5A{I) + SA{J) - 2SC{I, J)] 

(72) 
with the definition S{I) = Y^k^i^^^^^^)- I* is first 
obvious that A = is a iV times degenerate solution 
(we can arbitrarily fix the (5A(/)'s). If A 7^ 0, the first 
two equations define 5A{I) and 5C{I, J) as functions of 
dD(I,J), so that the third equation can be rewritten 
solely with the Z?'s : 

2n 2D 

iX-^){X+'-^)SD{I,J) = 

-(A 



2D 

— )iS{I) + S{J)). (73) 



It is now clear that either X — —^, which is satisfied for 
any SD{I, J) so that it is N{N — l)/2 times degenerate, 
or A ^ ~1^' ^^ which case one can further simplify the 
problem in : (A - ^)SD{I,J) = D{S{I) + S{J))/N. 
A = ^ is a solution, provided S{I) = for any / (this 
is true because A^ > 3). The restrictions imposed on the 
SD{I,J) imply the degeneracy is N{N — 3)/2. Let us 
now suppose A is not equal to any of the eigenvalues we 
already found. We have to solve an equation of the type 
f{X)SD{I,J) = giX){S{I) + 5(7)), so that the S(/)'s 
completely determine the SD(I, J). Here / and g are very 
simple functions of A (/(A) = A - 2D/N, and g = D/N 
does not even depend on A) , but in the general case with 
spin, one is also confronted to such a problem, that is why 
we keep quite general arguments. Summing this equation 
first over J, then over /, and noting S = ^^ S{I), we get 
two compatibility conditions : 

[/(A) -{N- 2)giX)]S{I) = g{X)S, and (74) 
[f{X)-2{N-l)giX)]S = 0. (75) 

These equations, in the spinless case read (A — D)S{I) — 
DS/N and (A - 2D)S = 0. Thus, either X = D with the 
condition S = 0, that is the degeneracy is (A^ — 1), corre- 
sponding to the free choice of (A'' — 1) variables S'(/), or 
X ^ D which automatically implies A = 2D (otherwise 
we get a zero eigenvector), and the degeneracy is 1 be- 
cause once S is given, everything is known. As D = 1, 



we thus have found (A'^— 1) eigenvalues equal to 1, and 1 
eigenvalue equal to 2. This is in agreement with the re- 
sults previously obtained in Sec. IV. Notice that we find 
exactly the eigenvalue 1, whereas in Sec. [V, we found 
(A^ — 1) eigenvalues 1 — Ae/N'^. This does not invalidate 
those previous results, because they were obtained in a 
perturbation expansion in e, which is in fact an expan- 
sion in £"i/A^"2 with 712 > ni. Thus our exact result 
simply shows that the (ni = l,n2 — 2) term we calcu- 
lated will be compensated by the {rii — 2,n2 = 2) term, 
once e is set equal to 1. Let us also notice that all the 
other eigenvalues we calculated exactly are either 0, or 
tend to in the A" — > 00 limit, as we expect from the 
RPA results, where only eigenvalues equal to 0, 1 and 2 
appeared. We furthermore recover the degeneracies pre- 
dicted in the RPA calculation, once the infinite A^ limit 
has been taken. 

We list all the eigenvalues and their d egen eracies we 
found for all six fixed directions, in T able |VII| . In agree- 
ment with the discussion in Sec. VC , the SU(A^) (-1-,+) 
direction is very similar to the spinless case with SU(2Af) 
symmetry. It is clear that both situations lead to the 
same eigenvalues. The degeneracies are not identical 
since in the spinful case, we always enforce the SU(2) 
spin symmetry of the pair interaction, which has no 
counterpart for the SU(2A^) spinless problem. It should 
be noticed that in the SU(A^), (-|-,-) case, we do not 
give exact results, but a 1/A^ expansion of the eigen- 
values. We only give the lowest order term in this ex- 
pansion, except for the eigenvalue approaching 1 in the 
infinite A^ limit, because this eigenvalue A ~ 1 + 2/N'^ 
is greater than 1 for any finite A'^. This explains the in- 
stability of the spin density wave fixed direction, that 
one observes in the numerical simulation. Furthermore, 
the absence of eigenvalue greater or equal to 1 in the 
SU(2), (-1-,-) case, shows that this fixed direction is sta- 
ble, which is confirmed by the numerical results, since 
after leaving the spin density wave regime, the flow ends 
in the "hot spots" regime. Let us flnally notice that the 
A = 1 -I- 2/N'^ + 0{1/N'^) eigenvalue is also compatible 
with the RPA result A^^^ ~ 1 - Ae/N'^, in the sense 
the 1/A^ term is zero in both results. The difference in 
the 1/A^^ terms comes, as in the spinless case, from the 
(s/N)'^ term that we have not computed. 



B. SIMPLE INTERPRETATION OF THE A = 1 
EIGENVALUES 

We have already interpreted the occurence of A = 1 
(or A = 1 + 0{1/N)) thanks to the results of Sec. |^. We 
will now show that these eigenvalues simply correspond 
to a lowering of symmetry, from SU(Af) to SU(Af — 1). 
In other words, it simply corresponds to the possibility 
of "turning off" one chain among A^. 

To see this, it is possible to compute the eigenvectors 
associated to these eigenvalues, and see their influence 
on the fixed direction. This can be checked to confirm 



TABLE VII: List of the eigenvalues for the six fixed directions. We give them in the form (eigenvalue, degeneracy), 
that all the results are exact, except for the SU(A'^), (+,-) case, where we give a 1/A^ expansion of the results. 
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Notice 



SU(iV) 
spinless 


SU(Ar) 


SU(iV) 
(+,-) 


SU(2) 
spinless 


SU(2) 
(+,+) 


SU(2) 
(+,-) 


{-2/N,N{N - l)/2) 


(-l/iV,iV(iV- 1)) 


(-3/iV, iV(iV - l)/2) 


(-1,1) 


(-1/2,2) 


(-1/2,2) 


(0,iV) 


(0,7V) 


{-V3/N,N) 


(0,2) 


(0,2) 


(0,3) 


(2/7V,7V(iV-3)/2) 


(l/7V,7V(iV-l)) 


(-l/7V,7V(7V-l)/2) 




(1/2,2) 


(1/2,2) 






(0,iV) 












(1/iV, iV(iV - l)/2) 












(V3/N, N) 












(3/iV, N(N - l)/2) 








(l,iV-l) 


(1,7V-1) 


(l + 2/7V^7V-l) 




(1,1) 




(2,1) 


(2,1) 


(2,1) 


(2,1) 


(2,1) 


(2,1) 



our statement. We here give a much simpler proof, re- 
stricting ourselves to the spinless case. Let us suppose we 
isolate the chain whose number is Iq. We assume a per- 
mutation symmetry among the N — 1 other chains. We 
still denote hy A ^ A{I), C = C{I, J) and D = D{I, J) 
the couplings involving only the A^ — 1 equivalent chains 
(/, J ^ Iq), and wc now denote a = A{Io), c ~ C(/o, J) 
and d ~ D{Io, J). Eq. ( [l7| ) then takes the form : 



NdtA 

Ndta 

NdtC 

Ndtc 

NdtD 

Ndtd 



(N - 2)D^ + d^ 
{N - 1)^2 

~d^ 

{N - 3)D^ + d^ + 2{A ~ C)D 

{N - 2)dD + d{a + A~ 2c) 



(76) 



These equations admit a few fixed directions. Among 
them, we of course have the SU(iV) fixed direction, with 
a^ A= 1-1 /N, and so on. We have the SU(7V-1) fixed 
direction with a = c = d — and A = 1 — 1/{N — 1), etc, 
for which the chain Iq is switched off. We in fact have 
one more SU(A^— 1) fixed direction, with A, C, D taking 
their SU(iV) values 1 - 1/N-l/N and 1, and a = A, 
c = C and d = —D. This comes from the symmetry of 
Eq. ^w under d -^ —d. It is also reminiscent of the RPA 
discussion of Sec. |lV| , where the components (f>{I) were 
only fixed in absolute value, their signs being free. 

If we study the linear stability of the SU(A'') fixed direc- 
tion, we get the following six eigenvalues : —2/N twice, 
twice, 1 once and 2 once. These eigenvalues are in 
perfect agreement wit h what we found in the previous 
section (see Table VII). In the {A, a, C, c, D, d) basis, the 
A = 1 eigenvector reads : 



{vA,Va,Vc,Vc,VD,Vd) 

2 2{N-1) 4 

iV' 



(77) 



N 



' N{N -2)' 



2 



N -2 



,1 



It is straightforward to check that va — vc + vd and 
Va 7^ Vc + Vd- Thus, if we add the vector defining 
the SU(A'') fixed direction and ev, with arbitrary e, the 
SU(A^— 1) symmetry will be maintained while the SU(Af) 
symmetry will be broken. Of course we cannot obtain 



the vector defining the SU(A^ — 1) fixed direction sim- 
ply by adjusting e, because the passage from SU(Af) to 
SU(A^ — 1) symmetry is non-perturbative. We can how- 
ever tell V does not connect the SU(iV) fixed direction to 
the second SU(7V - 1) one (with a = A, c = C , d = -D), 
because from the coefficients oiv, we see that a change in 
d implies a change in a of the same order of magnitude. 
Therefore, we are led to conclude that a perturbation 
along a A = 1 eigenvector should lower the symmetry, 
in a case where this eigenvector is marginally relevant. 
In order to confirm this picture, we have to go beyond 
linear perturbation theory which only tells us the A = 1 
eigenvectors are marginal. 



NON-LINEAR STABILITY ANALYSIS 



It is not possible to simply use Eq. (76) in order to per- 
form a non-linear stability analysis. Indeed, as our prob- 
lem is non-linear, the effect of a perturbation along two 
A = 1 eigenvectors (that lower the symmetry) cannot be 
determined by linear superposition. But neither are we 
forced to use the full RG Eq. ( p7| ) and perturb the fixed 
direction by all {N — 1) marginal eigenvectors. From the 
symmetry between the chains, it suffices to restrict the 
perturbation along two eigenvectors. Let us demonstrate 
this with the general notations introduced in Sec. IV. We 



will assume that we only perturb the fixed direction in the 
marginal subspace £i. From the flow equation, the cou- 
pling vector will also aquire components along the other 
eigenvectors, but these will be supposed to be negligible, 
which is reasonable since they will go to zero on a time 
scale that is short compared to the slow evolution in £i. 
We thus suppose gi{t) = f{t){ui + Vi{t)), with V € £i. 
We now plug this in the flow equation dtQi = Aijj-gjgu 
and get : 



dsV,^Vi[A,jkVj{s)Vk{s)] 



(78) 



where s{t) — \a[tc/{tc — t)] is the time scale adapted to 
the slow motion in £i, and where Vi is the projector on 
£i. It should be noticed that if the RHS vanishes, our 
approach is not valid anymore, and one has to take the in- 
fluence of the terms we neglected into account. If we now 



29 



TABLE VIII: a and 6 coefficients in the spinless and spin 1/2, 
(+,+) cases. 





SU(iV) 
spinless 


SU(iV) 
( + , + ) 


a 


4 


4 


{N-2)(N+2) 


{N-1)(N+1) 


b 


8 
(]V-2)^(]V+2) 




(N-2)(N-1)(N+1) 



decompose Vi{s) = X^QMal'S)^, , where a runs from 1 
to N — 1, and where the v's are marginal eigenvectors 
that lower the symmetry from SIJ{N) to SU(A^ — 1), we 
get : 

5.A*„ = J2 <M^(^) + E SfVMs)A*7(s) • (79) 

In this equation, the coefficients A and B are de- 
fined by VAAjkvfK[^^] = Eo ^^{"\ and V,[{A,,k + 

^*fcj>f^4''^] = Ea-Sf'^^l"^- We know from the Ap- 
pendix that RG flows cannot break a symmetry. This 
imposes severe conditions on A and B. If one perturbs 
the fixed direction along only one eigenvector [v'^^^), that 
is one singles out only one chain, the remaining SU(A^— 1) 
symmetry must be conserved, thus A^ — aSp^a- cl does 
not depend on a because of the permutation symmetry 
between the chains. If one now perturbs the fixed di- 
rection along two eigenvectors [v^^^ and «'■'''■'), and with 
the same perturbation strength along each vector, the re- 
maining SU(A^ — 2) symmetry as well as the permutation 
symmetry among the two vectors are conserved, so that 
-Bf ^ = b{5j3^a + ^7, a)- Again, b does not depend on a be- 
cause of the permutation symmetry between the chains. 
Eq. ([79|) now takes the much simpler form : 



dsfJ-a = a^^ 






(80) 



We have thus demonstrated that the non-linear evolution 
near the SU(A^) fixed direction is captured by only two 
numbers, a and b, which can most easily be determined 
by restricting oneself to two marginal eigenvectors, that 
is generalizing Eq. ([76|), with two singularized chains /i 
and l2- This tremendously reduces the dimension of the 
space we have to work with, since it is 11 in the spinless 
case, and 22 in the spin 1/2 case. The results are given 
in Table VIII. We see that there is no big difference 



between the spinless and spinning case. What is impor- 
tant is that in both cases, a is negative and behaves as 
1/N'^ when TV is large, and b is positive and behaves as 
1/N^. Let us see what happens for a perturbation along 
two eigenvectors. Eq. ( |80| ) gives : dsiJ,i = 6/ii(/i2 ^PMi) 
and dsH2 = ^M2(/^i — P/^2), with p = —a/b {p > 1 for 
A'' > 4). We have performed a numerical simulation 
of this flow (for the spinless values of a and b, and for 
A^ = 8), and the result is shown on Fig. ^. The dashed 
lines (with no arrows) are the lines of equation fj,2 = pfJ-i 
and fi2 = fJ-i/v- The first of these lines determines the 
set of points where the RHS of the flow equation for /ii 



/^2 n = 




FIG. 26: Structure of the flow for a perturbation along two 
marginal eigenvectors. The arrows indicate the direction of 
the flow along the trajectories. The dashed lines (with no 
arrows) are the lines of equation /i2 — p/ii and fi2 = Mi/p- 



vanishes. It thus delimitates two areas, one of growing 
fii (under this line), and one of decreasing /ii (above this 
line) . The same analysis applies for the second line (sym- 
metric of the first one, with respect to the first bisector), 
with /ii replaced by fi2- The arrows give the direction of 
the flow on the trajectories. It is important to notice that 
the point (^1 = 0,^2 = 0), corresponding to the SU(A^) 
symmetry, is a fixed point of the flow. The main con- 
sequence of our analysis is to show that this fixed point 
is unstable, except for the case where both fii and ^2 
are positive. If one had taken a perturbation along all 
marginal eigenvectors, the stability condition would have 
been that all A^ — 1 ^'s should be positive. The relative 
size of the stable subset thus goes to zero as 1/2^, so that 
in the large A^ limit, nearly all initial conditions lead to 
an instability of the SU(A'^) fixed direction. 

At this point, the reader would be right to point the 

) fixed 



IIC,thattheSU(A) (- 



fact that we said in Sec 
direction was numerically stable, and that our analytical 
results seem to refute the numerics. In fact, this fixed 
direction is found to be unstable, both analytically and 
numerically. At the beginning of this article, we didn't 
want to go into too much detail, so we did say the fixed 
direction was numerically stable. Anyway, this was not 
so much of a "lie" . Indeed, before analytical results were 
available, we believed it was stable, because the numer- 
ical simulations didn't show any sign of instability, for 
RG times of the same order than the ones needed to 
see the (+,-) instability. But once the analytical study 
pointed to instability, we tried to work with much big- 
ger RG times, and we indeed observed an instability. We 
should also mention that this instability could only be 
seen for quite small values of A^, because otherwise the 
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precision of the machine was not sufficient to detect it. 
Finally, let us mention the following point. In the case of 
the (+,-) fixed direction, the linear instability results in 
a direct transition between the SU(7V) and SU(2) fixed 
directions. In the two other cases, the non-linear insta- 
bility leads to a gradual elimination of pairs of chains, 
beginning by the centre (resp. the edges) for the spin- 
less (resp. (-f,-l-)) case. We shall not go further in the 
analysis of these transitions, and turn to our last section. 



VII. SYMMETRY RESTORATION AND THE 
LARGE N LIMIT 



The subject of symmetry restoration through RG flows 
is an interesting f|iipS|l,ir)p mWihich several authors have al- 



also find three fixed directions given by 



B6ff 



ready addressed .E3E3E3il3il3 We shall here take advan- 
tage of what we have said before, and show some impli- 
cations on symmetry restoration, in the large N limit. By 
symmetry restoration, we mean that the vector of cou- 
plings flows towards a symmetric flxed direction. We will 
focus on the case of the ajfiisotropic Gross-Neveu model, 
studied by Azaria et alel and which is related to the 
physics of coupled chains at half flUing. By contrast to 
the statements made by these authors, we believe that a 
careful analysis of one-loop RG flows is sufficient to indi- 
cate whether symmetry restoration is to be expected in 
the large N limit. 

Let us first set up our notations for this peculiar model 
(the general notations are the same as in the previous 
sections). We consider a ID Gross-Neveu model, with 
N types of fermions, numbered from to A'^ — 1. The 
0-fermions are singled out in order to get an anisotropic 
Gross-Neveu model. We thus get three possible density- 
density interaction terms. Denoting the size of the sys- 
tem by L, we get the following Hamiltonian : 



H = dx< i^ol''idf,ipo + 'ipal^idf.ipa 



+ f (^OV^O)' + ^iMo){^a^a) + jji^l^a^af \ , (81) 



with a = I, . . . ,N — 1 {N > 2), and where the sum 
over a is to be understood. The 1/A^ factors have been 
included to yield a good thermodynamical limit. Setting 
g — g/{2Tr) for the three coupling-Cjonstants, we get the 
following one-loop RG equations :C3 



' dtg = 
dtgc = 
. dtgo = 



-4 (1 -i)9'~ W 
-2 (2 -f). 95c -^5offc 



(82) 



(w, Uc, uo) 



N 



-1,0,0), 



(83) 



A{N -2) 

N 
(»,^e,^o)- ^^^_^^ (-l,±2,-l). (84) 

We shall only be interested in the last one, with Uc — 
2u — 2uq, corresponding to symmetric coupling con- 
stants involving all fermions. In the first case, ■00 
fermions do not interact. We have performed a linear 
stability analysis around the symmetric fixed direction, 
exactly as is Sec. VIA, and found the following eigenval- 



where t = Ih^Aq/v), as in Sec. Ill A. These equations 
admit a line of fixed points, given by g = g^ = 0. We 



ues for the matrix B : 

1 1 1 - 2/N 1 , , 

Notice that for A^ = 4, one gets the results already found 
by Lin.o We thus see that the symmetric fixed direction 
is locally stable, but it becomes marginal in the infinite 
N limit, as the A ~ 1 — 1/iV eigenvalue tends to 1 in this 
limit. This picture is the same as the one we developed 
in the previous sections, except it provides an example 
where an eigenvalue is strictly smaller than 1 when N is 
finite, and goes to 1 when N —f cx3. 

We shall consider the RG flows for two initial condi- 
tions. Let us first consider an initial condition, for which 
all couplings are negative and g = 2gc = go (we have 
not set hc{0) = h{0) because in this case, he — h a.t 
the beginning of the flow, and this makes it difficult to 
read the flgures). Notice that the sign of gc will only 
distinguish between the two fixed directions of Eq. (|84|), 
so that we can choose gc < 0. From the signs of the 
couplings and the RG flow equations, it is clear that the 
norm will quickly grow. The numerical flow in direction 
space has been computed with a large value of iV (10'*), 
and the results are shown on Fig. p^. Notice it is possi- 
ble to use such large values of N because here N is only 
a parameter, that does not determine the dimension of 
the coupling space, which is always equal to 3. One thus 
sees that in the very end of the flow (Fig. ^(a)), the 
symmetry is restored, since g = go and gc = 2g. How- 
ever, the physics will be given by the beginning of the 
flow. It is indeed clear from Fig. p^(c) that the norm 
will have exploded before the symmetry is restored. In 
Fig. |2^(6), the couphngs seem to be on a fixed direction. 
But we have seen there are only three fixed directions 
(Eqs. ( p3[ ) and (Q)). In fact, what we see here is a RPA 
like flow. Indeed, if one takes the infinite N limit in 
Eq. (JS^), one easily sees that we get an infinity of fixed 
directions, which will all be marginal. This is also con- 
firmed by the presence of the A ~ 1 — 1/iV — > 1 eigenvalue. 
This first example is very close, in spirit, to what we have 
dealt with in the previous sections of this article. 

Let us now consider another initial condition for which 
the flow passes in a close neighborhood of the line of 
fixed points. This is what Azaria et al have consid- 
ered. We thus choose g > and go < 0. The con- 
dition g > implies that the absolute values of both 
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FIG. 27; RG flow of difi'erent quantities, (a) : flow of the 
three couplings h, he, and ho represented respectively with 
squares, circles and triangles. (6) is a zoom on short time s 
of (a) . (c) shows the evolution of the norm, and (d) gives the 
relation between the two RG times s and t. 



g and gc will decrease in the early stages of the flow, 
and go < implies the instability of the intermediate 
fixed point. Again, the sign of gc is chosen to be neg- 
ative. Precisely, in the space of normalized couplings, 
we set : /i(0) = —2hc{Q) = — /io(0) (again we have not 
set /ic(0) = /i(0) to make the figures clear). The corre- 
sponding flows are shown on Fig. Eq. One thus sees that 
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FIG. 28: RG flow of different quantities, (a) : flow of the 
three couplings h, he, and /lo represented respectively with 
squares, circles and triangles. (&) is a zoom on short time a 
of (a) . (c) shows the evolution of the norm, and (d) gives the 
relation between the two RG times s and t. 

in the very end of the flow (Fig. p8|(a)), the symmetry 
is again restored. However, one has to look carefully at 
what happens at the beginning of the flow. It is clear 
from Fig. pq(6) and (c) that the flow gets close to the 
line of fixed points. It spends a very long time there 
(s < 50000). In the same time, t grows linearly with s 
(Fig. |28|(d)), with a proportionality coefhcient given by 
the i nver se of the norm of the couplings (remember from 
Sec. |IIB| that As = M{t)dd,). This shows that the sym- 
metry restoration will only take place at a very small 
energy scale (which goes to zero as exp(— TV)). Further- 
more, we have seen that the symmetric fixed direction 



is locally stable, but that it becomes marginal in the 
N -^ oo limit. Thus, the time needed to get close to 
this fixed direction, with a given accuracy, grows with 
N . But when the couplings are near this fixed direction, 
the norm grows nearly as fast as if the couplings were 
on the fixed direction, that is it explodes very quickly. 
We thus conclude that to observe a symmetry restora- 
tion (with a given accuracy) in the domain of validity 
of the one-loop RG, the couplings must be sent to zero 
as N goes to infinity. In other words, if the initial cou- 
plings are finite, and N is taken to infinity, the validity of 
the RG equations will break down before the symmetry 
is restored. We thus have found two indicators, given by 
the RG itself, showing that the symmetry restoration will 
not physically occur. 

For both initial conditions we have chosen, we have 
seen that the RG restores the symmetry. But this hap- 
pens in a regime where the pertubation theory is not 
valid anymore (and in the very far infra red for the sec- 
ond initial condition). What is important is that the RG 
itself tells us not to believe in the symmetry restoration. 



VIII. CONCLUSION 

In this article, we have revisited the old problem of 
interacting electrons with a short range repulsion in the 
weak coupling limit, but with a perfectly nested Fermi 
surface in two dimensions. Whereas a mean- field analy- 
sis predicts a spin density wave instability at low enough 
temperature, we have focused our study on the compe- 
tition between this instability and the onset of super- 
conductivity, which is known as the "parquet" problem. 
In purely one-dimensional sytems, this competition is a 
key ingredient in stabilizing the famous Luttinger liquid 
phase. In two dimensions, a careful analysist^ has shown 
that after a high-energy transient regime, the dominant 
effective couplings indeed correspond to a spin density 
wave state. According to this important work, the main 
effect of the competition between various instabilities is 
merely to lower the value of the critical temperature, in 
comparison to the mean-field estimate. In the present 
paper, extendipg our previous numerical study of RG 
flow equationSjE^ we have examined in greater detail the 
evolution of the effective couplins in the vicinity of the 
critical energy scale associated to the onset of long-range 
order. By contrast to the results of Zheleznyak et al,ta 
we find a third regime, in which the spin density wave 
instability slowly gives way to a d-wave superconductor. 
This final stage of the flow can be described as a "hot- 
spot" regime since it involves large couplings only when 
initial and final particles are located near the end points 
of the flat regions on the Fermi surface. 

We have analyzed the residual perturbation from 
particle-particle scattering on the mean-field like spin 
density wave regime, and we have found that its main ef- 
fect is indeed to gradually reduce the width of the Fermi 
surface regions with large 2kp couplings, as the typical 
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energy scale is reduced. However the same considera- 
tions show this effect is strongly dependent on the num- 
ber N of Fermi surface patches which are always intro- 
duced in the numerical solution of the RG flows. We 
have shown that the flow (in the space of directions for 
the coupling vector) slows down tremendously in the spin 
density wave regime, as N goes to infinity. Therefore, the 
magnitude of the coupling vector increases by a factor of 
order exp(7A'^^) (where 7 > 0) during this phase before 
the system reaches the third regime. So we conclude that 
initial couplings have to be extremely small (of the order 
of exp(— 7iV^)) if the cross-over from spin density wave 
to d-wave superconductivity is to be predicted in a reli- 
able way from the one-loop approximation. For real 2D 
systems with flat Fermi surface, N is in lactiittfinite, so 
the more traditional picture is recovered .r ^r W 1 

As a final note, we would like to make a brief remark 
on the various numerical computations of RG flow for 
the 2D Hubbck«(ir«i©del which have appeared over the 
past few yearsil'aBQU Since they introduce by numerical 
necessity a finite set of patches, do they suffer from the 
same limitation which has been discussed here ? We are 
inclined to say no in fact, since in most of these works, 
the bare couplings are not chosen to be very weak. The 
instability occurs then relatively quickly, so that all the 
subtle effects connected with the finiteness of N do not 
have time to influence the flow in a dramatic way. 
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imation). We have an important simple consequence of 
this which is : "if the effective couplings are invariant 
under the action of a generator of C at some scale, they 
will remain invariant at any other scale" . This combined 
property of RG flo ws and symmetry operations has been 
used in Sec. VIC, For any subalgebra ^ of £, we con- 



sider the subspace A'^ of all the two-particle interaction 
terms which commute with any element of A. We know 
that the RG flow leaves A'^ globally invariant in the sense 
that it takes any element of A'^ at scale Ag into another 
element of A"^ at scale A. Conversely, for any subspace 
S in the space of all two-particle interaction terms, we 
may define the subalgebra S"^ of all generators of £ which 
commute with any element of S. Clearly S is included 
in (S"=)^ We don't expect S = {S")" to hold for any 
subspace S, but if there is a subalgebra A oi C such that 
S = A", then S = (5"=)^ Qualitatively, the RG flow has a 
tendency to gradually eliminate all the non-relevant cou- 
plings with respect to the low-energy fixed point. Here 
a fixed point may be defined only in the projective cou- 
pling space, since the magnitude of the couplings is found 
to diverge at a finite energy scale, within the one-loop 
approximation. So the flow is confined to the vicinity of 
smaller and smaller subspaces S'„ as more and more irrel- 
evant directions disappear. This decreasing sequence of 
subspaces corresponds to an increasing sequence of sub- 
algebras A'^ ■ So the mechanism of symmetry generation 
at low energies seems to fit this simple-minded picture. 

Let us now turn to a more precise discussion of these 
points. Let us assume that our fermions can be in 
N' = 2N internal states (combining chain and spin de- 
grees of freedom), which corresponding creation opera- 
tors will be denoted by cj^ ^{k) and c^ ^{k). Let us con- 
sider the right moving branch, with the kinetic Hamilto- 
nian Hq = up Z^fc.al^-^p) : c\i^.a{k)cK,a{k) : . From this 
expression, it is clear that the following operators : 



APPENDIX : ONE-LOOP RENORMALIZATION 
GROUP AND CONTINUOUS SYMMETRIES 

As noticed before by Lin, Balcnts and Fisher jl^ the 
non-interacting Hamiltonian Hq for a system of coupled 
chains possesses a very large symmetry. For a system of 
N species of spin 1/2 fermions with the same Fermi veloc- 
ity, it is easy to exhibit a Lie algebra C of generators com- 
muting with Hq and isomorphic to so(4A^)pj^ so(4Ar)L. 
This very large symmetry is of course dramatically bro- 
ken by the interaction part of the Hamiltonian. But many 
people have noticed a general trend towards enhanced 
symmetries of the effective. .gnu tjI ings , as the RG flow 
runs towards low energies .E3Ejoc3'l3 This phenomenon 
is certainly related to an important property of RG flows : 
"The action of the generators of £ commutes with the RG 
flow" . If we denote the bare coupling vector at scale Aq 
by go, and the corresponding effective coupling at scale A 
by ^(go, Ag, A), then for any transformation a belonging 
to/:, we have: 0(3(50, Aq, A)) = ^(0(59), Aq, A). A proof 
of this statement is given below (for the one- loop approx- 



^R,a,/3 = E/c CR,a(^F + ^)4,/3(^F - k) 



D 



R,Q,/3 



E. 



-R,a 



(fc)cR,/3(A:) 



(86) 



commute with Hq. It is also well establishedlia that they 
form a closed Lie algebra isomorphic to C Eq D^^a.a ® 
so{2N'). A given Hermitian generator G in this Lie al- 
gebra operates on the fermion operators according to : 



[G,4^(fc)] 
[G,CR,^(fc)] 



A^,c- 



a,/3''R,a(^) + 
C^,/34,J2fcF 



B^,clj2kF 



a, (3^ 
k) 



k) 



DlAjk) 



(87) 



-A% and C% = -Bl^. 



From G = Gt , we get D^^^ 

Furthermore, the canonical anticommutation relations 
for the electron fields are preserved under the action of 
an infinitesimal transformation, which implies : 



B% + B^„, = Q 



^a,f3 



c^:« = o 
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Of course, we can construct similar symmetry generators 
for the left moving branch. In the presence of interac- 
tions, the huge £ — so{2N')yi ® so(2iV')L symmetry al- 
gebra of Hq is no longer preserved in general, but for 
some special choices of the pair interaction, some non- 
trivial subalgebras of £ may still commute with the full 
Hamiltonian. It is therefore useful to study the action 
of the generators of C on the pair interaction. In gen- 



eral, we may write the full bare Hamiltonian a,s : H — 
Ho + i j:k,k',, Eo,/3,,,5 ^(7, 5\a, f])cl^{k + q)cls{k' - 
q)cL,f3{k')cYi.a{k), where the coefficients F{j,6\a,P) 
parametrize the interaction. Taking G in £ as be- 
fore, and computing [G, H] yields the transformation law 
F ^ F + SF of the interaction function under the oper- 
ations of £. A simple computation gives : 



J 



6F{j, S\a, (3) = A«;^i^(A, 6\a, (3) + A^^Fi^f, \\a, f3) + dI^F{j, 5\\, p) + D^p^^Fi^i, 6\a, A) 



(89) 



r 



Note that in general some terms which do not conserve 
the total particle number such as c^c'c^c or c'ccc are also 
generated. We shall restrict ourselves to situations where 
those terms vanish. This corresponds to the following 
constraints : 



Bf^F{X,6\a,(3) 
Bl^F{j,X\a,P) 
CJ^^Fij,S\a,\)- 
[Cl,Fij,S\X,P) 



-Bl^FiX,6\j,P)^0 
B^.^F{j,X\a,S) = 
CJ,F{^,P\a,X) = 



(90) 



Now, the one loop flow Eqs. mdi) have the following form : 

dtF{-i,S\a,(3) ^ F^F 

= i^(7,/i|A,/3)F(A,5|a,/i) (91) 

~FiJ,S\X,^l)FiX,^l\a,P). 

We are now ready to make the two following state- 
ments : 



i) Suppose the action of a generator G of £ on an in- 
teraction function F satisfies the conditions of Eqs. (pC|). 
Then F * F also satisfies these conditions. 

ii) Choosing G and F as in i), the action of G com- 
mutes with an infinitesimal RG transformation, namely : 
{F + 6F) * {F + 6F) = F * F + S{F * F) + 0{SF^). Or 
equivalently : S{F =*= F) = F ^ 5F + 5F '^ F (we use the 
notation SF as in Eq. (p9) above) . It is a simple exercise 
to check that i) and Mjhold. We just need the last two 
relations of Eq. ( ^8|) for i) and the first one for ii). 

The consequences of i) and ii) are immediately 
extended to finite (instead of infinitesimal) RG trans- 
formations. So the property that a generator G in £ 
preserves the total particle number conservation is valid 
at any stage along a RG flow trajectory, provided it 
holds for the initial couplings. Similarly, invariance 
of an interaction function along a transformation G 
propagates well along any RG flow trajectory. 
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